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We study the CMB observables in axion monodromy inflation. These well-motivated scenar- 
ios for inflation in string theory have monomial potentials over super-Planckian field ranges, 
with superimposed sinusoidal modulations from instanton effects. Such periodic modulations 
of the potential can drive resonant enhancements of the correlation functions of cosmological 
perturbations, with characteristic modulations of the amplitude as a function of wavenumber. 
We give an analytical result for the scalar power spectrum in this class of models, and we 
determine the limits that present data places on the amplitude and frequency of modulations. 
Then, incorporating an improved understanding of the realization of axion monodromy infla- 
tion in string theory, we perform a careful study of microphysical constraints in this scenario. 
We find that detectable modulations of the scalar power spectrum are commonplace in well- 
controlled examples, while resonant contributions to the bispectrum are undetectable in some 
classes of examples and detectable in others. We conclude that resonant contributions to the 
spectrum and bispectrum are a characteristic signature of axion monodromy inflation that, 
in favorable cases, could be detected in near-future experiments. 
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1 Introduction 



Inflation [H El [3] is a successful paradigm for describing the early universe, but it is sensitive 
to the physics of the ultraviolet completion of gravity. This motivates pursuing realizations 
of inflation in string theory, a candidate theory of quantum gravity. Considerable progress 
has been made on this problem in recent years, so much so that the most pressing task, 
particularly in view of upcoming CMB experiments, is to learn how to distinguish various 
incarnations of inflation in string theory from each other and from related models constructed 
directly in quantum field theory. 

Fortunately, the additional constraints inherent in realizing inflation in an ultraviolet- 
complete framework can leave imprints in the low-energy Lagrangian, and hence ultimately in 
the cosmological observables. In favorable cases, a given class of models may make distinctive 
predictions for a variety of correlated observables, allowing one to exclude this class of models 
given adequate data. 

One decisive observable for probing inflation is the tensor-to-scalar ratio, r. A promising 
class of string inflation models producing a detectable tensor signature are those involving 
monodromy [I], in which the potential energy is not periodic under transport around an 
angular direction in the configuration space. The first examples [I] involved monodromy 
under transport of a wrapped D-brane in a nilmanifold, and a subsequent class of examples 
invoked monodromy in the direction of a closed string axion [5]. 

The axion monodromy inflation scenario of [5] is falsifiable on the basis of its tensor signa- 
ture, r m 0.07. However, primordial tensor perturbations have not been detected at present, 
while the temperature anisotropies arising from scalar perturbations have been mapped in 
great detail [6]. One could therefore hope to constrain axion monodromy inflation more ef- 
fectively by understanding the signatures that it produces in the scalar power spectrum and 
bispectrum. Characterizing these signatures is the subject of the present paper. 

As we shall explain, the potential in axion monodromy inflation is approximately linear, 
but periodically modulated: each circuit of the loop in configuration space can provide a 
bump on top of the otherwise linear potential. Modulations of the inflaton potential with 
suitable frequency and amplitude can yield two striking signatures: periodic undulations in 
the spectrum of the scalar perturbations, and resonant enhancement [7] of the bispectrum. 
Let us stress that the presence of some degree of modulations of the potential is automatic, 
and is an example of the situation described above in which traces of ultraviolet physics 
remain in the low-energy Lagrangian. We do not introduce modulations in order to make 
the scalar perturbations more interesting. However, it is important to examine the typical 
amplitude and frequency of modulations in models that are under good microphysical control, 
in order to ascertain whether well-motivated models produce signatures that can be detected 
in practice. 

To achieve this, we first investigate in detail the realization of axion monodromy inflation 
in string theory. We compute the axion decay constants in terms of compactification data, 
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we assess the importance of higher-derivative terms, and we estimate the amplitude of mod- 
ulations for the case of Euclidean Dl-brane contributions to the Kahler potential. We also 
identify a potentially-important contribution to the inflaton potential, arising from backre- 
action in the compact space, and we present a model-building solution that suppresses this 
contribution. 

We find that detectable modulations of the scalar power spectrum and bispectrum are 
possible in models that are consistent with all current data and that are under good mi- 
crophysical control. In fact, we find substantial parameter ranges that are excluded not by 
microphysics, but by observational constraints on modulations of the scalar power spectrum. 

The organization of this paper is as follows. We begin in £j2] by describing the classical 
evolution of the homogeneous background in axion monodromy inflation with a modulated 
linear potential. We then solve, in £j3j the Mukhanov- Sasaki equation governing the evolution 
of scalar perturbations, giving an analytical result for the spectrum in terms of the frequency 
and amplitude of the modulations of the potential. Next, we briefly discuss the bispectrum 
and express the amplitude of the non-Gaussianity in terms of the model parameters. We then 
present, in £|4| an analysis of the constraints imposed on axion monodromy inflation by the 
WMAP5 data (for prior work constraining similar oscillatory power spectra, see e.g. [8l f9| fTOl. 
ITlj H2| \13\ EH])- Then, in ^5] and £j6j we present a comprehensive analysis of the constraints 
imposed by the requirements of computability and of microphysical consistency, including 
validity of the string loop and a 1 perturbation expansions, successful moduli stabilization, 
and bounds on higher-derivative terms. In £j7]we combine the observational and theoretical 
constraints, with results presented in figure [7| 

1.1 Review of axion monodromy inflation 

In this section we will briefly review the motivation for axion monodromy inflation, as well as 
the most salient phenomenological features. We will postpone until §[5] a more comprehensive 
discussion of the realization of this model in string theory. 

Inflation is sensitive to Planck-scale physics: contributions to the effective action arising 
from integrating out degrees of freedom with masses as large as the Planck scale play a 
critical role in determining the background evolution, and hence the observable spectrum of 
perturbations (see [15] for a review of this issue). A central problem in inflationary model- 
building is establishing knowledge of Planck-suppressed terms in the effective action with 
accuracy sufficient for making predictions. The most elegant solution to this problem is to 
provide a symmetry that forbids such Planck-suppressed contributions. Because invoking 
such a symmetry amounts to forbidding couplings of the inflaton to Planck-scale degrees of 
freedom, it is important to understand this issue in an ultraviolet-complete theory, such as 
string theory. 

One promising mechanism for inflation in string theory involves the shift symmetry of 
an axion. Axions are numerous in string compactifications and generally enjoy continuous 
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shift symmetries a — > a + constant that are valid to all orders in perturbation theory, but 
are broken by nonperturbative effects to discrete shifts a — > a + 1. As noted in [5J, the 
shift symmetries of axions descending from two-forms are also broken by suitable space-filling 
fivebranes (D5-branes or NS5-branes) wrapping two-cycles in the compact space. 

In axion monodromy inflation [5], an NS5-brane wrapped on a two-cycle S breaks the 
shift symmetry of the Ramond-Ramond two-form potential C*2, inducing a potential that is 
asymptotically linear in the corresponding canonically normalized field 0, 

v = i?<j>, (l.i) 

with n a constant mass scale. Inflation begins with a large expectation value for the inflaton, 
oc J s Ci 3> 1, and proceeds as this expectation value diminishes; note that the NS5-brane, 
like any D-branes that may be present in the compactification, remains fixed in place during 
inflation. As argued in [5], this gives rise to a natural model of inflation, with the residual 
shift symmetry of the axion protecting the potential from problematic corrections that are 
endemic in string inflation scenarios. 

In this paper we perform a careful analysis of the consequences of nonperturbative ef- 
fects for the axion monodromy scenario. Such effects are generically present: specifically, 
Euclidean D-branes make periodic contributions to the potential in most realizations of axion 
monodromy inflation. However, the size of these contributions is model-dependent. It was 
shown in j5] that there exist classes of examples in which nonperturbative effects are practi- 
cally negligible, but we expect - as explained in detail in §6.5| - that in generic configurations, 
periodic terms in the potential make small, but not necessarily negligible, contributions to 
the slow roll parameters. 

Therefore, it is of interest to understand the consequences of small periodic modulations 
of the inflaton potential in axion monodromy inflation. In this paper we address this question 
in two ways: first, in §[2}§|4j by studying a phenomenological potential that captures the 
essential effects; and second, in £|5]and £[6j by investigating the ranges of the phenomenological 
parameters that satisfy all known microphysical consistency requirements dictated by the 
structure of string compactifications in which axion monodromy inflation can be realized. 



2 Background Evolution 

In this section we will study the background evolution of the inflaton in the presence of small 
periodic modulations of the potential. We will focus on modulations in axion monodromy 
inflation with a linear potential, but our derivations are easily modified to account for other 
models with a modulated potential. We will denote the size of the modulation by A 4 , and 
write our potential as in [3], 



V (0) = /i 3 + A 4 cos 



/ 



+ bf cos 



/ 



(2.1) 



4 



where we defined the parameter b 



-4ij. The equation of motion for the inflaton is then 



+ 3H<j) + /x 3 - /i 3 6 sin 



0. 



(2.2) 



To solve (2.2), we begin with two approximations. Monotonicity of the potential require^] 
b < 1, and as we will see in £|4| for the case b < 1 observational constraints in fact imply b ^ 1. 
This suggests treating the oscillatory term in the potential as a perturbation. Furthermore, 
the COBE normalization implies that <fi 3> M p during the era when the modes that are 
observable in the cosmic microwave background exit the horizon. This allows us to drop 
terms of higher order in M p /<p. 

Under these conditions, it is straightforward to solve for the evolution of the homogeneous 
background. Expanding the field as <fi = O + b(j)\ + 0(b 2 ), the equations of motion of zeroth 
and first order in b become 



(2.3) 



+ V3/2 



3 • I r0 
ji° sin ( — 



(2.4) 



where we have neglected terms of higher order in M v j and we have made use of the slow 
roll approximation for O Using equation (2.3), we can rewrite equation (2.4) with O as an 
independent variable instead of t, yielding 



3 
2 



?o sin 



(2.5) 



where primes denote derivatives with respect to 0o- F° r the period of interest, in which the 
modes now visible in the CMB exit the horizon, it is a good approximation to neglect the 
motion of 0o everywhere except in the driving term. The inhomogeneous solution is then 
given by 



Mt) = f 



(2 + 3/ 



2\2 



36/ 2 



-(2 + 3/ 2 )sin 



Mt) 

f 



+ 6/0* cos 



Mt) 

f 



(2.6) 



1 The case of non-monotonic potentials may also be interesting. On the one hand, for sufficiently large 
b > 1, it may be possible to realize chain inflation [TSJ [T71 [TS] in our model. In this scenario, the inflaton 
would tunnel from minimum to minimum, with the universe expanding by less than one third of an e-fold 
per tunneling event. This requires a more careful analysis, and we will leave this for future studies. On the 
other hand, for b > 1 the model essentially turns into a small-field model of inflation because the inflaton gets 
trapped at the peaks for a large number of e-folds. It seems hard to distinguish this from other models of 
small field inflation, but it may be interesting to take a closer look at this as well. 

2 In approximating sin(<^//) ~ sin(0 o //) on the right hand side of (2.4 1, we have assumed not only that 
b <C 1 but also that b<j>i/f -C 1. As we will see from the solution (2.7), <j)\ is of order / 2 0*. Hence the mild 
assumption <C 1 justifies this approximation. 
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where 0* denotes the value of the field O & t the time at which the pivot scale k* exits the 
horizon. Assuming 60 e-foldings of inflation, this happens around 0* ~ 11M P . For decay- 
constants / obeying / > M p /10, there is less than one oscillation in the range of modes that 
are observable in the cosmic microwave background, leading to an uninteresting modulation 
with very long wavelength. We will thus make the additional assumption that / <C M p . 
Assuming that 0o ^ M p and /< 1, using the slow roll approximation for 4>o(t), and working 
to first order in b, the solution thus becomes 



M) = Mt) + bMt) = Mt) + 6 / 1+ 3 ( 3^)2 

with (f>o(t) given by 



— sin 



Mt) = 



Mt) 
f 



2/3 



+ 3/0* cos 



Mt) 

f 



(2.7) 



li s/2 (t-U) 



(2.8) 



In the absence of oscillations, i.e. for 6 = 0, axion monodromy provides a model of large 
field inflation that is easily studied using the slow roll expansion. Assuming for concreteness 
that the CMB scales left the horizon 60 e-foldings before the end of inflation, we are interested 
in the perturbations around 0* ~ 11M P . After imposing the COBE normalization, one finds 
that CMB perturbations are produced at a scale V 1/4 ~ 7 • 10~ 3 M P ~ 1.7 • 10 16 GeV with 
a spectral tilt n s ~ 0.975 and a tensor-to-scalar ratio r ~ 0.07. For reference, the Hubble 
constant during inflation is then H ~ 2.8 • 10~ 5 M P ~ 6.8 • 10 13 GeV. 

One can then ask what happens once the oscillations are switched on, i.e. when b ^ 0. 
It turns out that the effect on the number of e-foldings is negligible as long as 6 <C 1. Hence 
the inflationary scale is well-approximated by the slow roll analysis. On the other hand, the 
detailed properties of the perturbations are very different from the slow roll case and cannot 
be calculated in that expansion. We turn to this issue in the next section. 



3 Spectrum of Scalar Perturbations 



Having understood the background evolution, we are now in a position to calculate the power 
spectrum in axion monodromy inflation. One might be tempted to do this by brute-force 
numerical calculation, but we find it more instructive to have an analytic result. We will 
show that under the same assumptions made in calculating the background evolution, i.e. 
slow roll for <j>o(t), 0o ^> M p , f M p , and to first order in 6, the scalar power spectrum is of 
the form 



A 2 



TZ 



(k) 



A 2 



n 



1 + Sn s cos 



f 



A 2 



-!+r 



n(fc/fc») 



(3.1) 



where the quantity A^(A;*) parameterizes the strength of the scalar perturbations and will be 
introduced in detail in the next subsection. The second equality is valid as long as 5n s <^ 1, 
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and 5n s is given by 



5n, 



12b 



where 



y/{l + WW) V 8 



fc = a/02 - 2 In fc/fc, ~ 0* 



7T / 7T 

COth I ~ ~ ; ) /0h, 



2/0, 



In kj k* 



(3.2) 



(3.3) 



is the value of the scalar field at the time when the mode with comoving momentum k exits 
the horizon. 



In £3.1 we will give a derivation of this result that makes no further approximations. In 



£3.2 and £3.3 we will present two additional derivations of (3.1) that are valid only as long as 



<C 1 but that lead to a better understanding of the relevant physical effects behind the 



power spectrum (3.1). Let us at this point briefly summarize the scales that will be relevant 



for our discussion in the next subsections. 



Given the potential (2.1 ), the time frequency of the oscillations of the inflaton is lu = (j)/f. 



This is also the time frequency of the oscillations of the background. Perturbations around 
this background can be quantized in terms of the solutions of the Mukhanov- Sasaki equation, 
assuming an asymptotic Bunch-Davies vacuum. Every perturbation mode with comoving 
momentum k oscillates with a time frequency k/a that is redshifted by the expansion of the 
universe until the mode exits the horizon and freezes when k = aH. 

Then, if H < uj < M p , every mode will at a certain time resonate with the background, 
as stressed by Chen, Easther, and Lim in [7]. Using the slow roll equation of motion and the 
COBE normalization, 



3H<j> ~ -V'(<t>) , 



~-eV, F~5-10~ 7 eM 



p ■ 



the requirement H < uo < M p can be re-expressed as 



to 
H 

UJ 

Ml 



M 2 
p 



> 1 



2eV 1 



< 1 



(3.4) 

(3.5) 
(3.6) 



hence defining a range of values for the axion decay constant / for which resonances occur. 



Using V2e ~ M p /0, ~ .09, we obtain 2.4 • KT 6 < jL- < 0.09. We will show in ^5 



and £6 



that 



/ falls in this range in a class of microphysically well-controlled examples. 

Going beyond our approximations, the model also predicts a small amount of running of 
the scalar spectral index, of order 10~ 4 , from terms of higher order in the M p /0 expansion. 
Furthermore, 5n s develops a very mild momentum dependence. We will neglect these effects 
because these will most likely not be observable in current or near-future CMB experiments. 
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3.1 Analytic solution of the Mukhanov- Sasaki equation 

We begin our study of the spectrum by choosing a gauge such that the scalar field is unper- 
turbed, 50(x, t) = 0, and the scalar perturbations in the spatial part of the metric take the 
form 

Sgytx, t) = 2a(t) 2 ft(x, t)S lj . (3.7) 

The quantity TZ(pc, t) is a gauge-invariant quantity and in the case of single-field inflation is 
conserved outside the horizon. It is closely related to the scalar curvature of the spatial slices, 
but we will not need its precise geometric interpretation at this point. 

The translational invariance of the background and thus the equations of motion govern- 
ing the time evolution of the perturbations make it convenient to look for solutions of the 
linearized Einstein equations in Fourier space. One defines 

/d 3 k 
[n k (ty^a(k) + n k (ty e -^a*(k)] , (3.8) 

where k is the comoving momentum, and k is its magnitude. The rotational invariance 
of the background ensures that H k (t) can depend only on the magnitude of the comoving 
momentum but not on its direction. Directional dependence can only be contained in the 
stochastic parameter a(k) that parameterizes the initial conditions and is normalized so that 

(a(k)a*(k')> = 5(k - k') , (3.9) 

where the average denotes the average over all possible histories. With this ansatz, the 
Einstein equations turn into an ordinary differential equation, the Mukhanov- Sasaki equation, 
governing the time evolution of 1Z k (t). We will use it in the forr^ 



d 2 TZ k 2(1 + 2e + 5) dK 



+TZ k = 0, (3.10) 

dx z x ax 

where x = — kr, with the conformal time r given as usual by r = J* Outside the horizon, 
i.e. for x <C 1 or equivalently k/a <C H, the quantity 7Z k (x) approaches a constant which 
we denote by 7Z k °^ . In terms of 7Z k °^ we define the primordial power spectrum for the scalar 
modes as 

2 -'f. (3.1D 

To evaluate this quantity, it will again turn out to be sufficient to solve to first order in b. We 
therefore expand the slow roll parameters, 

e = e + e 1 + O(b 2 ), (3.12) 



3 We use the same definitions for the slow roll parameters as in [19] , i.e. e = — -§2, S = - . S is related to 
the Hubble slow-roll parameters r\ = e/eH by S = rj/2 — e. The other slow-roll parameters that are sometimes 
used are ey = (V /V) 2 /2 and rjv = V" /V. When the slow roll expansion is valid they arc related to the 
Hubble slow-roll parameters by ey = e and rjv — 4e — r/. 
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5 = 5 Q + 6 1 + 0(b 2 ) . 



(3.13) 



For the background solution (2.7), the first-order terms are given by 



3bf 



0*[1 + (3/0*) 2 ] 



cos | yj + (3/0*) sin 



(3.14) 



36 



[i + (m*) 2 

We now consider an ansatz of the form 



sin 



/ 



(3/0*) cos 



/ 



(3.15) 



n k = n 



k.o 



Here the index uq on the Hankel function, Hol(x), is given by uq = | + 2eo + <5o, (?(x) is a 
perturbation of order 6, and 7^°q is the value of 7Zk(t) outside the horizon in the absence of 
modulations, i.e. for 6 = 0. To be explicit, it is given byQ 



(3.16) 



K 



(o) 
k,0 



Ml 1 

6 k 3 / 2 



(3.17) 



where 0^ ~ 0* — ln ^ //fc * once again is the value of the scalar field at the time the mode with 
comoving momentum k exits the horizon. The quantity of interest to first order in b is then 



K 



(o) 



K 



(o) 
k.O 



l + 2Rec/(0) 



K 



,2 Re 9 (0) 



K 



(o) 
A,0 



2Rcg(0) 
^ \ ln(fc/fc«) 



(3.18) 



Our ansatz automatically solves the equation of order b°. To first order in b and in the slow 
roll parameters, the Mukhanov-Sasaki equation leads to an equation for g(x) of the form 



#9 
dx 2 



2dg^ 
x dx 



+ g = 2e ix (2e 1 + S t 



(3.19) 



In writing this equation, we have dropped terms of order O(beo, b5o), which amounts to setting 
Vq = 3/2. Next, we notice that e\ is suppressed relative to 8% by a factor Since we are 
interested in the regime 4- -C 1, we can thus drop the term proportional to e\ on the right 



hand side of equation (3.19). Furthermore, it turns out to be convenient to rewrite Si using 



trigonometric identities. Ignoring an unimportant phase, one finds 



Si 



3b 



V 1 + (3/0*) 2 



cos 



/ 



(3.20) 



4 As mentioned earlier, we will ignore the running of the scalar spectral index, but it may be worth pointing 
out that the information about the running is contained in this formula. 
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It will be convenient to write 4>q{x) as 4>o(x) 
r(x) = Re(g(x)), equation (3.19) becomes 



d 2 r 2 dr 



66 



ln(fc/fc») , lnx 



lnx 

COS(X) cos I — + — — 

/ J<p* 



k + ^p. Introducing 



(3.21) 



cte 2 x y/i + (3/0*) 2 

The solution to this equation can be found e.g. using Green's functions. We are particularly 
interested in the inhomogeneous solution at late times, i.e. in the limit of vanishing x. Using 
more trigonometric identities, we find that the solution in this limit can be brought into the 
form 



r(0) 



cos I -f + P(ffa) 
l + (3/0*) 2 \f 



where /3(/0*) is an unimportant phase that we will ignore, and X is the integral 



dxJi(x) J_i (x) x«> 

2 2 



(3.22) 



(3.23) 



Written in this form, the integral can be recognized as a Weber-Schafheitlin integral and can 
be done analytically (see e.g. [20J). One finds 



\I\ 



' 7T / 7T 

coth 



2/0* 



Combining equations (3.18), (3.22) and (3.24), we finally obtain an expression for Sn 



5n x 



2r(0) 



126 



cos 



V 7 ! + (3/0*] 



coth 



7T 



2/0* 



(3.24) 



(3.25) 



Once again, this derivation is valid to first order in b and assumes slow roll for (j)o(t), 0o ^ M p , 
and / M p . In particular, it makes no use of an /0* 1 expansion, although this 
approximation will be needed in the derivations in £372] and §3.3[ A comparison between our 
analytical result for Sn s as a function of for a fixed value of b and the result of a numerical 
calculation using a slight modification of the code described in [2T] is shown in Figure [l] 



3.2 Saddle-point approximation 

As we have seen in the last subsection, it is possible to calculate the power spectrum ana- 
lytically to first order in b, assuming slow roll for 4>o(t), O ^> M p , and / M p , but the 
derivation sheds little light on the physics behind the results. To get a better understanding, 



it is instructive to look at the integral (3.23 ) more explicitly. For this purpose, it is convenient 



to separate X into its real and imaginary parts, X 



X c + iX S} with 



X r 



J, 



, (smx — xcosx) cosx / In x 
dx cos 



dx 



x* 

(sinx — x cosx) cos a; flnx 
„ sin 



;r- 



f4>* 



(3.26) 
(3.27) 
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Figure 1: The solid line is the analytical result for 5n s as a function of /, for b = 0.08, while 
the dots are the numerical result obtained from an adaptation of the code used in 



For ranges of the axion decay constant such that /</>* <C 1, these integrals can be done in 
a stationary phase approximation. Using trigonometric identities to rewrite the products 
of trigonometric functions appearing in the integrands into sums of trigonometric functions 
with combined arguments, one finds that the stationary phase occurs at x = jhr- Expanding 
around the stationary point and performing the integral as usual, one finds to leading order 



in 



X. 



X, 



7T 



sin 



7T 



COS 



1 + ln(2/^) n 

f<t>* 4_ 

l + ln(2/^) rr 

fd>* 4 



which leads to 



IXI 



1r ~\~ 1 « 



7T 



(3.28) 
(3.29) 

(3.30) 



This agrees with our previous result, equation (3.24), as long as /</>* <C 1. We have not 



only reproduced our earlier results, however: we also learn that at least for small 



2kf<p* 



the 



Up to the factor of two in 



integral is dominated by a period of time around f - 
the denominator, this corresponds to the period when the frequency of the oscillations of 
the scalar field background equals the frequency of the oscillations of a mode with comoving 
momentum The stationary phase approximation thus captures a resonance between the 
oscillations of the background and the oscillations of the fluctuations, and is good as long as 
fcf)* C 1, i.e. as long as the resonance occurs while the mode is still well inside the horizon. 



°This factor of two can be understood from momentum conservation, as will become clear in £3.3 
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One might suspect that this has an interpretation in terms of particle production, and we 
shall make this more precise in what follows. 

Recall that our ansatz for IZk was given in (3.16 ), where g(x) is the solution of the equation 

d 2 g 2 dg 



dx 2 x dx 



+ g = 2e lx 8 1 



with 81 again given by 



36 

V 1 + (3/0*) 



cos 



t lax 

+ ~rr 



(3.31) 



(3.32) 



and initial conditions given by lim g(x) = and lim g'(x) = 0. As we have just learned, 

x— +00 x^oo 

the effect of the driving term can be ignored long after the resonance has occurred, i.e. for 
x <C 2/^~j^] This implies that at late times, g(x) must be a solution of the homogeneous 
equation which can be written as 



/ , / Zx^H^hx) I + c 



'3/2^ 



(3.33) 



where cjj. are momentum dependent coefficients. The solution for equation (3.31) can also 
be written explicitly as 



g(x) 



(x cosx — smx) 



2e iy (cosy + ysmy) 



Si 



(3.34) 



+ (cosx + xsinx) 



2e w (smy — ycosy) 



y 



Si. 



For x <C we can take the lower limit in the integrals to zero and this can be brought into 
the form 



g(x) = \{1 2 + il 1 ) (i^xlHi 1 ^ + \{Z 2 - zX 1 ) [~i^K 
where the integrals X 1 and X 2 are given by 



(2), 
3/2 < 



(3.35) 



X 1 



X 2 



66 



V 1 + (3/0*)' 
66 

V 1 + (3/0*)' 



e iy (cos y + y sin y) lnx 



cos 



/ 



+ 



e iy (sin y — y cos y) /0a,. In x 



cos 







/ 



+ 



(3.36) 



(3.37) 



6 One should note that this is not because the driving term goes to zero, but because its frequency becomes 
too high for the system to keep up with it. 
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In the saddle point approximation these evaluate to 



1 = iZ 



6b 



71" -(4>k l+ln2/0 t 7r 



(3.38) 



Combining equations (3.16), (3.35), and (3.38), we finally find that the curvature perturbation 
for x <C 



takes the form 



x) 



with ci "* given, up to an unimportant momentum-independent overall phase, by 



(3.39) 



,(-) 



6b 



Vi 1 + (3/0* 



7T 



(3.40) 



One might now interpret the coefficient of the negative frequency mode as a Bogoliubov 
coefficient that measures the amount of particles with comoving momentum k being produced 
while this mode is in resonance with the background. It seems hard to make this precise as 
one really is comparing mode solutions of different backgrounds rather than mode solutions 
of different asymptotically Minkowski regions in the same background. 



Equation (3.39) also shows that instead of starting in the Bunch-Davies state and then 



following the mode through the resonance, one may start the evolution after the resonance 
has occurred but use a state that is different from the Bunch-Davies state, which is similar 
to what is considered in [22], |23j |8j |9j [10]. The departure from the Bunch-Davies state is of 
course quantified by c^~\ 



3.3 Particle production and deviations from the Bunch-Davies state 

Here we will deal with a conceptual question that generically arises in inflationary models 
with oscillations in the scalar potential. Driven by the background motion of the inflaton, the 
oscillating contributions constitute a time-oscillating perturbation to the Hamiltonian of the 
system. Now, perturbations oscillating in time will generically induce transitions, in our case 
from the original vacuum state to some excited states. This implies that the vacuum state of 
the full system will deviate from the Bunch-Davies vacuum of the homogeneous background 
inflationary evolution. We will now estimate the resulting quantity of particle production 
and relate the result to the derivation of the scalar power spectrum given in the preceding 
sections. 

To lowest order the oscillating perturbation is given by 

AH^ = ^V"(Mt))-H\ (3.41) 
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implying that the lowest-order transitions will be from the vacuum |0) to two-particle states 
\p, —p). As the physical momentum p = k/a corresponding to a given comoving momentum 
k is exponentially decaying in the inflationary regime, any two-particle state with given co- 
moving momentum \k, —k) will be in resonance with the oscillating perturbation only for a 
short period of time which we will have to estimate in due course. 

In transforming the Hamiltonian of the fluctuations into Fourier space 

= \&%+ \p 2 HI+ \v"(Mt)) ■ 54, (3.42) 

we find that the system takes the form of a perturbed harmonic oscillator with eigenfrequency 
u p — p = \p\ for each momentum mode 8p separately, 

H[H P \ = \&% + + \v"{m) ■ Hi • (3.43) 



2 -rp 2 p -rp 2 

vr - 



A_ff(2)[50 p ] 

Now we compare this to the perturbed harmonic oscillator in one-dimensional quantum 
mechanics, 

H= l -x 2 + l -uolx 2 + ^6u{t) V . (3.44) 
In going to dimensionless variables q,p we can write this as 

H = \u (p 2 + g 2 ) + ^£g 2 , (3.45) 

AHP> 

where for our case of a periodic perturbation periodic with frequency u we have 

5uj(tf = 5uj 2 cos(cjt) . (3.46) 

We want to determine the time-dependent transition matrix element in time-dependent per- 
turbation theory for a periodic perturbation. To do so, we first write the perturbation in 
standard form for time-dependent perturbation theory as 

A#( 2 )(t) = ^l! g 2 = ^tfuiut + g-i-t) = p r e i*t + e -fa,tN ) (3 _ 47) 

in the notation of equations (40.1) through (40.9) of [21]. The Hamiltonian and the transition 
matrix elements can be written in terms of creation and annihilation operators a and using 
q = (a + a^)/v2 and p = —i(a — af)/y2. Then, canonical quantization of the unperturbed 
part yields a discrete spectrum \n) of eigenstates with energy spectrum E n = u>o(n + 1/2). 

If we compare this with our actual case above, we see that for each momentum mode 5(pp, 
q and p are replaced by appropriate dimensionless fields 5<pp and 11^. In complete analogy 
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to the simple quantum mechanical oscillator, there will be a tower of discrete states \n) p with 
energies E n>p = u p {n + 1/2) = p(n + 1/2). In particular, |2) p labels the two-particle state 

2uj. 

thus have for the perturbation in our actual case 



\p, —p) which has energy difference AE 2)P = 2u p = 2p with respect to the ground state. We 



AH {2 \t) 



5u(ty 2 

-oft 



2u r 



4uj, 



,,jt _|_ e ~Wt\ = p (Jut _|_ e ~wt\ 



(3.48) 



For the transition matrix element one then finds 



(p, -p\AH^\0) = F 20 (e^ + e- 
Here we have used that 



■iiot \ 



with F' 



6u' 



\p, -p) 



20 



,-h2 



a 2 ( a + a) 



|o) 



5u 2 



V2 



-|o> 



and 



{0\a 2 (a + (J) 2 \0) = (0\a?((f) 2 \0) 



(3.49) 



(3.50) 



(3.51) 



If the energy of the two-particle state E 2p = 2k /a were not too close to the perturbation 
frequency u, we could use time-dependent perturbation theory with the above matrix element 
and obtain the first order transition probability Po^2fc, 

2 



Po- 



•2k 



dt'{k,-k\AH {2 \t')\0)e 



(f) 2 +- 2 + (fr-- 2 



'2k\ 2 



2\F- 



cos(2ujt) 



20 | 



Sou 4 



.(f) 2 --. 

(f) 2 +- 2 +kf) 2 



cos(2a;t) 



16(fc/a) 



2Ai\ 2 



(f) 



(3.52) 



where cu 2 o = E 2)P =k/ a — E 0jP=k / a = 2k /a. This gives the resonance line feature characteristic 
of transition processes. 

However, as for any given k the physical momentum and frequency k/a will decrease 
extremely rapidly with 1/a, we can approximate the amount of transition happening in the 
short time interval At res during which the two-particle state of given k is in near-resonance 
uj m 2k /a. Close to resonance, time-dependent perturbation theory breaks down (visible in 
the singularity of the above result for uj = 2k /a); however, for periodic perturbations one can 
solve the Schrodinger equation of the coupled two-state system exactly |24j. One finds that 
on resonance the transition probability is 



Po^2k = \ [1 - cos (2ftt)} = l - 



cos 



2V2k/< 



-t 



where Q 



20 



(3.53) 
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That is, near resonance the system effectively oscillates with frequency 2Q = 2 



between the vacuum and the two-particle state. 

We now have to estimate the time At res during which a two-particle state of comoving 
momentum k stays in near-resonance. We will follow the analysis in [7] and look at the 
interference terms induced between the cos(o;t) perturbation and the exp^u^t) periodicity 
of the interaction matrix element in (3.52). We note that, on the one hand, the two-particle 
state with frequency 2k /a = uj — Auj stays in resonance with the perturbation with frequency 
u) only for a time roughly estimated to be (for the relative phase shifting from — 7r to 7r) 

Ati ~ — . 3.54 
Auj k j 

On the other hand, in the inflating universe it takes very roughly a time 

At 2 ~ ^ (3.55) 

to change the frequency of the two-particle state from, say, to + Auj to uj — Auj. Equating the 
two provides us with the effective duration of near-resonance, 



At res = Ati = At 2 ~ 2i/vr— if" 1 . (3.56) 

V uj 

Plugging this into the above transition result and remembering that near resonance k/a ~ 
uj/2, we get 



1 



cos 




(3.57) 



Now, in our case above we see that p = k/a in terms of comoving momenta k, and further 
5uj(t) 2 = V"(Mt)) = 4r cos (^^] = 5u 2 cos(ut) , u = ^- . (3.58) 



P \ f 

Noting that in our scenario of interest we have H < uj and that 5uj H, we can expand 
the argument of the cosine around zero. If we then plug in the microscopic definitions of the 
quantities 5u 2 = A 4 // 2 and to = H/(f(f), we get 

^tt A 8 ff_7rA7^_7r 9A 8 /^ _ 367r^ 

- 2 f4uj2H2 ■ ^ 2 /2#4 2 /2/i 6 2 g ^- 0J ^ 

where 0* ~ 11M P denotes the vev of the inflaton field around 60 e- foldings before the end of 
inflation. 

Next, because P^ 2 k characterizes the transition probability to the two-particle states, 
it may be related to the negative frequency Bogoliubov coefficient that relates the out- 
vacuum to the in-vacuum. Specifically, the out-vacuum is specified by the modes 

Ukiput) = c^'u-k + c^ + \k , (3.60) 
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whereas the original Bunch-Davies in-vacuum had modes 



Uk{in) = Uk , S + \in) = 1 



(3.61) 



We therefore find that 



(3.62) 



In comparing these results with the general treatment of the Mukhanov-Sasaki equation 



above, we see by looking at (3.39) and (3.18) that we can identify 



Uk 



U-k 



(3.63) 
(3.64) 



and thus from (3.18) we conclude that 

2Reg(x) 2Re(c< + >c<->) nl , 

- ~ 2 c ( " 



5n s 



cos 



x^O 



COS 



(3.65) 



which agrees with the general result (3.25) in the appropriate limit u > H and 5u <C H, 
corresponding to /0* < 1, where coth (7r/2/0*) — > 1. 

Note that in calculating the transition probability we lose information about the phase 



of the transition matrix element as given in (3.49). Therefore, if we estimate the population 



coefficient from \/Po->2fcj we get only an estimate for \c^ \ without the phase information. 
A more complete derivation using the full information in the transition matrix element should 
also yield the information about the phase as derived in the previous subsection. 

Thus, we see that in the regime of rapid oscillations, < 1, the induced 5n s is due to 
a time-localized deviation from the Bunch-Davies state, which may be interpreted as being 
due to resonant bursts of particle production happening well before a given mode leaves the 
horizon during inflation. 



3.4 Bispectrum of scalar perturbations 

We start by reviewing how resonance can drive the production of large non-Gaussianity during 
inflation, as proposed in [7] . We then present an estimate for the size of the non-Gaussianity 
for the model (2.1). 

The three-point function can be calculated as [25] 

(7e(r,k 1 )7e(r,k 2 )^(r,k 3 )) = -% r<[^(r,k 1 )^(r,k 2 )^(r J k 3 ), J Hj(T / )]>adT / , (3.66) 

J Tn 
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where Hj is the interacting part of the Hamiltonian. Hj was calculated for a generic potential 
(see e.g. [2SJ CTJ) at cubic order in the perturbations; it takes the form 



- / d 3 x 



ae 2 nn' 2 + ae 2 n{dny - 2eR!{&lZ){d X ) 



where d denote space derivatives, 



X = a 2 ed 2 1Z 



(3.67) 



(3.68) 



and we used the Hubble slow-roll parameter r\ = e/(eH) = 2(e + 5) because formulas in this 
subsection are simpler in terms of r] than in terms of 5. 



We would like to stress that (3.67) is exact for arbitrary values of the slow roll parameters 



e and rj. Substituting Hi into (3.66) produces six terms, plus an additional term coming from 



a field redefinition. For the modulated linear potential (2.1), e is small, as in standard slow 



roll inflation. On the other hand, contrary to the standard slow-roll approximation, fj can 
be much larger than e 2 . This suggests that the leading term comes from the ef) term in the 
Hamiltonian]^] Hence we have [HI [7] 



n 



Ui(T end ) X 



7~end / A 

dreri'a 2 \ uI{t)u* 2 {t) —u* 3 {t) + sym ) 5 3 (K)(2tt) 3 + c.c. . 



As in [7], we parameterize the non-Gaussianity as 



G(ki, k 2 , k 



s) x3 



{k\k 2 h) 



^(K)A^(2vr) 



(3.70) 



(3.71) 



where K = ki + k2 + k3. We take as an ansatz for the shape of the non-Gaussianity for our 
modulated linear potential 



Following [7J and comparing (3.70), (3.71) and (3.72), we obtain the estimate 



Sr 



3 771 



7 In (3.9) of this term was written as 

12 



8H, 



2 d (/) 



p 2 dt \ 24>p vy 

which can be reduced to the term in (3.671 using H' = —<f> 2 /2. 



(3.72) 



(3.73) 



(3.69) 
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where we have again used the notation r] = 7] + br]i + . . . . Using the background solution 
obtained in £l2| it is straightforward to find 



/i 3 6b 



cos (Jj j + (3/0*) sin (^j^j 



(3.74) 



3^/[l + (3M) 2 ] 

It is not hard to convince oneself that in the region of parameter space where f res > 1 and 



6 < 1, the second term in (3.74) is always negligible, i.e. 3/0 <C 1. Hence our estimate for 
the non-Gaussianity is 

/»=^=?»^r. (3-75) 



4(/0) 3 / 2 4 \HJ 

where we remind the reader that u = <p/f. As we will often refer to this equation, let us 
pause and comment on it. The resonant non-Gaussianity vanishes when the modulation is 
switched off, i.e. for 6 = 0. It is inversely proportional to some power of / (depending on 
which quantity is held fixed). Hence the smaller the axion decay constant /, the larger the 
non-Gaussianity. On the other hand, as we will see in §[5j there are theoretical lower (as well 
as upper) bounds on /, so that the non-Gaussian signal cannot be made arbitrarily large. 

No complete analysis of the observational constraints on resonant non-Gaussianity has 
been performed to date (however, see [26]), and such an analysis is beyond the scope of 
the present work. Based on a rough comparison with known shapes of non-Gaussianity, 
we estimate that f res > 200 might be at the borderline of being excluded by the current 
data, while f res < 1 would be difficult to detect in the next generation of experiments. A 
comprehensive analysis of the detectability of resonant non-Gaussianity is a very interesting 
topic for future research. 



4 Observational Constraints 

In the last section, we derived the theoretical predictions of axion monodromy inflation for 
the primordial power spectrum. We will now use these predictions to compare the model with 
the five-year WMAP data [Bj. While the data in principle allows for a variety of statistics 
to be extracted, we will limit ourselves to the most fundamental one, the angular power 
spectrum. The reason for this is that the data is not now adequate for the polarization data 
or the three-point correlations to place meaningful additional constraints on the model. This 
will change as soon as the Planck data becomes available, and will be an interesting problem 
especially given the unusual shape of the non-Gaussianities the model predicts. 

For the benefit of the less cosmologically-inclined reader, we now briefly summarize the 
basic observables relevant to our analysis. In the ideal scenario, in which a full-sky map is 
available, the temperature of the cosmic microwave background as a function of the position 
in the sky can be expanded in spherical harmonics as 

T (™) = ^2 a emY em {n) (4.1) 
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The theoretical counterparts of these measured expansion coefficients, which we will denote 
a^, should be thought of as random variables satisfying a (possibly only nearly) Gaussian 
distribution. Each realization of these coefficients corresponds to a possible history of the 
universe. In the Gaussian case, all the information about the theory is contained in the two- 
point correlations of these, as the odd n-point functions vanish, and the even n-point functions 
are sums of products of the two-point functions. Assuming an isotropic background, the two- 
point correlations must take the form 



a 



tm a t'm') - C t baib mm i , (4.2) 



where the brackets denote an average over all possible histories or equivalently (by the ergodic 
theorem) all possible positions. The a^'s themselves, being random variables encoding initial 
conditions, cannot be predicted from a given cosmological model, and only the multipole 
coefficients, Ce, encoding their correlations are of interest. These multipole coefficients Ci 
can be estimated from the measured expansion coefficients a^ m via 

c ? y = ^TiJ2M 2 ■ ( 4 - 3 ) 

m 

For noiseless, full-sky CMB data, these provide an unbiased estimate of the true power spec- 
trum in the sense that the average of the analogously defined quantity for the a^'s satisfies 

<^ sky ' th > = 2*TT /£ |afj2 ) = Ce ' (44) 

Since there are only 2£ + 1 modes per £, even for the ideal noiseless full-sky map the estimate 
of the multipole coefficient has the cosmic variance uncertainty 

(4.5) 



' ,osky,th ^ 
Kj c — L>f 



Ci J 2£+l 

In a more realistic setting with noise and sky cuts, this estimator is no longer unbiased 
and more sophisticated estimators have to be used. The current state of the art is to use a 
pixel-based maximum likelihood estimator for low I (specifically, for I < 32), and a pseudo-C^ 
estimator for higher I. For details we refer the reader to [27J and references therein. 

After this quick review of the basic relevant quantities, let us describe our analysis. We 
work on a grid of model parameters. For each point on the grid, we compute the theoretical 
angular power spectrum with the publicly-available CAMB code [2"8|l29j^j using the primordial 
power spectrum derived in the previous section in the form 



Ai(k) = AUK) ( - ) . (4.6) 



8 Of course, we modify the CAMB code to calculate all the multipole coefficients rather than calculating 
some and interpolating. 
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The likelihood for a given theoretical power spectrum is calculated with a modified version 
of the WMAP five-year likelihood code that is now available on the LAMBDA webpage [50] . 
The power spectrum in our model contains additional parameters beyond those of the WMAP 
five-year ACDM fit (namely, {Qsh 2 , Q c h 2 , Q\, r, n s , A 2 ^} and the marginalization parameter 
{Asz})- The additional parameters are 8n s , f and a phase A<p. This phase parameterizes both 
our uncertainty in the number of e-folds needed, which originates in our poor understanding 
of reheating, and a microscopically determined phase offset in the sinusoidal modulation of 
the scalar potential arising in the string theory construction. 

We fix the value of the scalar spectral index n s = 0.975. As in any model of large-field 
inflation, the spectral index is a prediction of the model that depends only on the physics of 
reheating and, correspondingly, on the total amount of inflation since the observable modes 
exited the horizon. The value we choose corresponds to the situation in which the pivot 
scale exits the horizon 60 e-folds before the end of inflation. The results turn out to be fairly 
independent of the precise value chosen for the scalar spectral index and we could have chosen 
the value corresponding to any number of e-folds between 50 and 60. We fix {Q c h 2 , r, A sz ) 
to the WMAP five-year best-fit values for the ACDM fit. We allow /, 5n SJ f2 B /i 2 , Aip to vary 
on the grid, and we also marginalize over the scalar amplitude { A^} in the likelihood code. To 
obtain Figure [2j we thus marginalize over {fi#/i 2 , A 2 ^} and over the unknown phase Ay?, while 
we fix {fl c h 2 ,fl\,r, Asz}, as we expect at most mild degeneracies between these parameters 
and the primordial ones. 

The grid consists of 16 equidistantly spaced points in fl B h 2 between fl B h 2 pa 0.0212 and 
Q B h 2 pa 0.0266, 128 equidistantly spaced points in 5n s between 5n s = and 5n s = 0.44, 
512 logarithmically spaced points in the axion decay constant / between / = 9 x 10~ 5 and 
/ = 10 _1 , as well as 32 points for the phase Atp between A(p = —n and A(p = n. This leads 
to a grid with a total of 33,554,432 points. The analysis was run on 64 of the compute nodes 
of the Ranger supercomputer at the Texas Advanced Computing Center. The compute nodes 
are SunBlade x6420 blades, and each of the nodes provides four AMD Opteron Quad-Core 
64-bit processors with a core frequency of 2.3 GHz. 

The resulting 68% and 95% contours in the 5n s — f plane are shown in the left plot of 
Figure [2j To convert the resulting observational constraints on 5n s as a function of / into 
constraints on the microscopic parameter bf as a function of /, we make use of equation (3.25 ). 
The resulting 68% and 95% contours in the bf-f plane are shown in the right plot of Figure [2] 
Roughly, the results can be summarized as bf < 10~ 4 for / < 0.01 at 95% confidence level. 
Our best fit point is at a rather small value of the axion decay constant, / = 6.67 x 10~ 4 , and 
a rather large amplitude for the oscillations, 5n s = 0.17. The fit improves by A% 2 ~ 11 over 
the fit in the absence of oscillations. The corresponding angular power spectrum is shown in 
Figure [3j 
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Figure 2: This plot shows the 68% and 95% likelihood contours in the Sn s -f, and 6/-log 10 / 
plane, respectively, from the five-year WMAP data on the temperature angular power spec- 
trum. 
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Figure 3: The left plot shows the angular power spectrum for the best fit point / = 6.67x 10 , 
and 8n s = 0.17. The right plot shows the angular power spectrum for the best fit point 
together with the unbinned WMAP five-year data. 



The improvement can be traced to a better fit to the data around the first peak. We 
would like to stress, however, that we do not take this as an indication of oscillations in 
the observed angular power spectrum. Similar spikes in the likelihood function occur quite 
generally when fitting an oscillatory model to toy data generated with the conventional power 
spectrum without any oscillations, because the oscillations fit some features in the noise. The 
polarization data could provide a cross check, but we find that it is presently not good enough 
to do so in a meaningful way. 
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Let us say a few words motivating the necessity of marginalizing over fig/i 2 and A(p. There 
is a known degeneracy in the angular power spectrum between VLsh 2 and n s , as changing Vtsh 2 
changes the ratio of the power in the first and second acoustic peaks, which to some extent 
can be undone by changing the spectral tilt n s . In our case we do not vary n s , but we add a 
sinusoidal contribution to the standard power spectrum. It is intuitively clear that by doing 
so we can change the ratio of power in the first and second acoustic peak by choosing the 
right oscillation frequency (controlled by /) and phase Ayj, leading to a degeneracy between 
Qsh 2 and Sn s at least for a certain range of /. 

The most straightforward way to demonstrate this degeneracy between f2#/i 2 and 5n s 
arising for certain 'resonant' values of / is to present a likelihood plot in the VLb-5u s plane for 
a value of / for which the degeneracy is clearly visible. An example is shown in the plot on 
the left side of Figure |4j It shows that marginalizing over VLbIi 2 is necessary to obtain correct 
exclusion contours on 5n s and /. That marginalization over the phase is necessary can easily 
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Figure 4: These plots show the 68% and 95% likelihood contours for the five-year WMAP 
data on the temperature angular power spectrum in the Sn s - Qsh 2 plane and Sn s - A0 plane 
for an axion decay constant of / = 3 x 10 -2 and / = 1.5 x 10~ 2 , respectively. 



be seen from a likelihood plot in the 8n s -A<p plane. This is shown in the plot on the right 
side of Figure [4} 

We have also performed a Markov chain Monte Carlo analysis for the model using the 
publicly available CosmoMC code [31], [32]. While the Monte Carlo has the advantage that 
it is less computationally intensive than a grid when varying all cosmological parameters, the 
likelihood function for oscillatory models turns out to be rather spiky, making the Monte 
Carlo hard to set up, because the chains tend to get trapped in the spikes. 

To some extent this can be overcome by taking out the problematic regions or increasing 
the temperature of the Monte Carlo. When run on parts of the parameter space where the 
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Monte Carlo runs reliably, we found agreement with the grid-based results shown above. 
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Figure 5: This figure shows a triangle plot for some of the parameters that were sampled in 
a Markov chain Monte Carlo for an axion decay constant of / = 10~ 2 . The contours again 
represent 68% and 95% confidence levels. 

The most problematic direction to sample is that of the axion decay constant, /. We show 
the result of one of our chains for / = 0.01 in Figure [5j The plot shows marginalized one- 
dimensional distributions and two-dimensional 68% and 95% confidence level limits for the 
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most important ACDM parameters as well as Sn s and A0. In the Monte Carlo, we sampled 
the parameters 8n s , A<f>, all parameters of the ACDM except the scalar spectral index, as well 
as the Sunyaev-Zel'dovich amplitude. 



5 Microphysics of Axion Monodromy Inflation 

In §1.1 we briefly reviewed the properties of axion monodromy inflation, focusing on the 
description in effective field theory. For a general characterization of the signatures of the 
scenario, the phenomenological model of §1.1 was sufficient. However, the phenomenological 
parameters /, //, A are in principle derivable from the data of a string compactification, and 
as such they obey nontrivial microscopic constraints: the ranges and correlations of these 
parameters are restricted by microphysics. 

We should therefore determine the values of the phenomenological parameters allowed in 
consistent, computable string compactifications. We will begin by reviewing the string theory 
origin of axion monodromy inflation, both to set notation and to highlight the properties 
most relevant in constraining the parameters /, /i, A. For concreteness we will restrict our 
attention to a specific realization of the scenario, in 03-07 orientifolds of type IIB string 
theory, with the Kahler moduli stabilized by nonperturbative effects. Our considerations 
could be generalized to other compactifications, but the numerical results would differ. 



5.1 Axions in string theory 

Let us first review the origin of the relevant axions. Our conventions and notation are sum- 
marized in appendix [Aj Consider type IIB string theory compactified on an orientifold of a 
Calabi-Yau threefold X. Let the forms u 1 be a basis of the cohomology H 2 (X, Z), normalized 
such that J 2 us 3 = 5j J \2n) 2 a' , where £/ are a basis of the dual homology H 2 (X, Z). The RR 
two-form C*2 gives rise to a four-dimensional axion via the ansatz^ 



C 2 = ^c I {x)u I , (5.1) 

where i is a four-dimensional spacetime coordinate. The ten-dimensional Einstein-frame 
action [33] that follows is 

/jio 9sy~9E I Jri ,2 / jio 9sy~9E nun i j a' jj' /<- \ 

X 2(2tt)V 4 1 21 = J X 12(2tt)V 4 ^ % c i d » c *>iPi<f9 E 9% ■ (5.2) 

Notice that the axions only have derivative couplings, and hence enjoy a continuous shift 



symmetry at the level of the classical action. In £5.2.3 we will recall the origin of this 



9 The factor of 2ir is introduced so that the four-dimensional axions c/ have periodicity 2ir, as can be seen 
via S-duality from the world-sheet coupling i J £?2/(27ra'). Notice that in our conventions C2 and uj 1 have the 
dimensions of length-squared, while cj are dimensionless. 
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symmetry and explain how it persists to all orders in perturbation theory and is broken 
by nonperturbative effects. 

Upon dimensional reduction, one finds a relation between the four-dimensional reduced 



Planck mass M p and a', 



a'M 2 p 



Ve 

7i 



(5.3) 



where Ve is the Einstein-frame (dimensionless) volume of the Calabi-Yau X measured in units 
of l s = 2^\fa' ; . The decay constant of the canonically normalized axion is then 



P 



9s 



M 2 

p 



f uj A *u> 



(5.4) 



48tt 2 V e [ (2vr) 6 a /3 

The present definition of the axion decay constant differs by a factor of 2tt from that in [5], 
i.e. f here = 2irf there . As a consequence our canonically normalized axion has periodicity 2irf, 
consistent with (2.1). 



5.2 Dimensional reduction and moduli stabilization 
5.2.1 Four-dimensional data of 03-07 orientifolds 

Now we consider how to stabilize the compactification in a setup that will allow inflation. We 
focus on the KKLT scenario for moduli stabilization [31]. We assume that the complex struc- 
ture moduli, the dilaton, and any open string moduli have been stabilized at a higher scale, 
and we concentrate on the remaining closed string moduli (specifically, the remaining moduli 
are those descending from hypermultiplets). The Af = 1 supersymmetric four- dimensional 
theory resulting from dimensional reduction of type IIB orientifolds was worked out in detail 
in [33]. We are interested in orientifold actions under which the holomorphic three- form Q 
of the Calabi-Yau manifold is odd, so that the fixed-point loci are 03-planes and 07-planes. 
The cohomology decomposes into eigenspaces of the orientifold action, 

#(*-.») = #;> s © H r _f • (5.5) 

We therefore divide the basis oj a, A = 1, . . . , ft, 1 ' 1 into u> a , a = 1, . . . h 1 ^ 1 and u> a , a = 1, . . . hi . 
Working out the sign of the orientifold action on the physical fields, one finds that the two- 
forms C*2 and B2 are odd, and should be expanded in terms of the u a . Grimm and Louis 
[35] have derived the Kahler coordinates on the corresponding moduli space, i.e. the proper 
complex combinations of fields that appear as the lowest components of chiral multipletsj 1 "] 

G ° - i^-'j) < 5 - 6 > 

T a = i Pa + ^c Q/3 yV + jC abc G b (G - Gf (5.7) 

10 We use the same notation as [35 with two exceptions: we rescale T a as T^ ere = (2/3)T^ here , and we add 
a factor of (27r) _1 in the definition of G a such that the fields c a and b a have periodicity 2ir. See appendix |A| 
for more details on our conventions. 
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where p a comes from the RR four-form C4 integrated over some orientifold-even four-cycle 



S a , with a 



/iV 1 ; and c a and b a come from the RR and NS-NS two-forms Co and B? 



integrated over some orientifold odd two-cycle S a with a 
potential is given by 



hi . The tree-level Kahler 
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K = log 



2\ogV E 



(5. 



where the (dimensionless) Einstein-frame volume Ve of the Calabi-Yau manifold is defined in 



(A. 4). The dependence of this Kahler potential on the multiplets (5.6) and (5.7) cannot be 



written down explicitly for a generic choice of the intersection numbers cjjk- The implicit 
dependence is given by writing the (Einstein-frame) volume in terms of two-cycle volumes v a 



K = log 



2 log 



6 



c aPy v a (T,G)v f> (T,G)v''(T,G) 



(5.9) 



Then one has to solve (5.7) for v a and substitute the result into the above Kahler potential. 



The Kahler potential is a function of Ve, and hence is a function of v a , and in turn of 
r a = Re T a and Im G, but does not depend on Re G and Im T a (as can be seen by taking the 



real part of (5.7)). One might be tempted to conclude that c enjoys a shift symmetry but 
that b does not, but, as we will explain in §5.2.3 , both fields have shift symmetries. 



The tree-level superpotential W does not depend on the multiplets (5.6) and (5.7). In fact 



it depends on the complex structure moduli and the dilaton, which we assume have already 
been stabilized by fluxes. Therefore we will take Wq to be a discretely tunable constant. 



5.2.2 Nonperturbative stabilization of the Kahler moduli 

Let us now proceed to consider nonperturbative effects. We follow the KKLT strategy [31] for 
the construction of a de Sitter vacuum. We assume that each four-cycle T a is wrapped either 
by a Euclidean D3-brane or by a stack of D7-branes giving rise to a four-dimensional gauge 
theory that undergoes gaugino condensation]^] This results in the following four- dimensional 
superpotential: 



W 



W + A a e~ aaTa 



(5.10) 



a=l 



where A a will be treated as constants, as they depend on the complex structure moduli, 
which we have assumed to be stabilized; a a = 2ir/N a , with N a the number of D7-branes in 



n We assume that the axio-dilaton t = Co + ze~^ is already stabilized by fluxes at r = i/g s and we write 
down the dilaton-dependent part of the Kahler potential only to keep track of factors of g s . 

12 In general, Euclidean D3-branes or D7-branes will wrap some linear combinations T a of the cycles ap- 



pearing in (5.7), rather than the basis cycles T a themselves, but for simplicity we will suppress this issue. 
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the stack; and N a = 1 for the case of a Euclidean D3-brane. We can find a supersymmetric 



minimum by solving for the vanishing of all the F-terms: for the h 1 ^ 1 even Kahler moduli via 

r d Ta V E 







D a W = d Ta W + Wd Ta K = -A a a a e~ aaTa - 2W- 



V E 



- W 



2V f 



and for the Zii' 1 odd moduli via 



D a W = d G aW + Wd G «K = -iW 



AttVe 



(5.11) 



(5.12) 



where in both cases in the last step we used the chain rule and the definitions of G a and T a 



in terms of two-cycle volumes v a . The condition (5.11 ) is simplified if we first solve for lmT a , 
which gives 



a a Im T a = 9 Aa - 8 Wo + k a ir , k a e Z . 
Then we are left with the set of real equations for each a, 

(±l) a \A a \a a e- a ^ = d Ta K ^\W \ + ^(±l) a |^| e -°^ j 



(5.13) 



(5.14) 



where (±1) depends on the value of k in (5.13). As long as the orientifold-even four-cycle 



Kahler moduli are defined as in (5.7), then d^ a K = —v a / (2Ve) < for every a. Now we prove 



that in (5.14) the minus sign has to be chosen for every a in order to have a supersymmetric 



solution. First we notice that the sign of the right hand side does not depend on a, so k ( 
and hence (±1) Q have to be the same for every a. If we choose the positive sign in (5.14) 



the quantity in brackets in the right hand side is manifestly positive. Then the two sides 
of the equation have opposite signs and no (compact) solution exists. To summarize, the 
minimization of Im T a boils down to taking all A a real and negative and Wq real and positive 
or the other way around,^ 

Concerning (5.12), an obvious solution is given by b a = for every a. As argued in [5], 



an inflationary model with a 6-type axion as the inflaton will generically suffer from an eta 
problem, and we will therefore focus on a c-type axion. 



5.2.3 Nonperturbative breaking of axionic shift symmetries 

Axionic shift symmetries are central to this paper, so we will now explain how they originate 
and how they are ultimately broken by nonperturbative effects. First, let us recall the classic 



13 



We notice that if one chooses as Kahler variable a linear combination of the T a defined in (5.7 1, as is done 



e.g. in the large volume scenario with Swiss-cheese Calabi-Yau manifolds, then the sign of dx a K can depend 
on a. In this case, the minimization of ImT Q boils down to taking Wq real and positive and A a real with 
Sign(A Q ) = Sign(9j' Q _R'), up to multiplying W by an overall phase. 
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result [3"6"l 137] establishing the shift symmetry to all orders in perturbation theory. Consider 
the axion b = J E 5/(27ra / ), where B is the NS-NS two-form potential and E is a two-cycle 
in the Calabi-Yau manifold. The vertex operator representing the coupling of b to the string 
worldsheet is [36] 

V(k) = ^ J <P^exp(ik'X^)y^d a X' t dpX v B ltv (X). (5.15) 

At zero momentum, this coupling is seen to be a total derivative in the worldsheet theory. 
Therefore, the axion b can only have derivative couplings (which vanish at zero momentum), 
to any order in sigma-model perturbation theory. Notice that the genus of the worldsheet did 
not enter in this argument, so the axion shift symmetry is also valid to all orders in string 
perturbation theory. 

This argument fails in the presence of worldsheet boundaries (i.e., D-branes), and also 
fails once worldsheet instantons, or D-brane instantons, are included. In axion monodromy 
inflation, both sorts of breaking play an important role, as we shall now explain. 

First, the introduction of an NS5-brane wrapping a curve S a creates a monodromy for 
the axion c a , spoiling its shift symmetry and inducing an asymptotically linear potential [5]. 
Specifically, the potential induced by the Born-Infeld action of the NS5-brane (obtained by 
S-dualizing the Born-Infeld action of a D5-brane) is 

V(c a )= \ /2 ^l* + (2ng s c a )^ (5.16) 

where iy/aS is the size of S a and e captures the possibility of suppression due to warping. For 
c a ^> 1, this potential is linear in c a , or in the corresponding canonically normalized field, 
which we denoted by <fi in the preceding sections. Let us remark that the square root form of 
the potential can be important at the end of inflation and also makes a small change in the 
number of e-foldings produced for given parameter values, so that in a model that includes 
a specific scenario for reheating, the square root structure should be incorporated as well. 
As we have not invoked a concrete reheating scenario, for our purposes the linear potential 
suffices, but one must still bear in mind that this form is not valid for small <fi. 

As we will explain in detail, the D-brane instantons involved in moduli stabilization in- 
troduce sinusoidal modulations to the linear potential. We will work exclusively in a regime 
in which the breaking by wrapped branes dominates over the nonperturbative breaking, al- 
though we remark in passing that the complementary regime might be interesting for realizing 
models involving repeated tunneling. 

The breaking of the b shift symmetry by Euclidean D-branes (or by gaugino condensation 
on D7-branes) is slightly subtle, so we will address it briefly. As we remarked above, b appears 
quadratically in the classical Kahler potential, which seems to contradict the statement that it 
enjoys a shift symmetry at the perturbative level in the absence of boundaries. However, there 
is no contradiction: the shift symmetry of b is true at constant two-cycle volumes v and not 
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at constant four-cycle volumes T. To see this, suppose that there is a single Kahler modulus 
T, so that the superpotential is of the form (5.10) with h 1 ^ 1 = 1. The Kahler potential is then 
[35] 

K= -3M p 2 log (T + T -db 2 ), (5.17) 

with d a constant. In the absence of a nonperturbative superpotential, a suitable simultaneous 
shift of T + T and b is a symmetry of the scalar potential of this system; under such a shift, 
the two-cycle volumes v are invariant. However, this symmetry is spoiled by the nonpertur- 
bative term in W, because the superpotential and the scalar potential are no longer invariant. 
Therefore, in a scenario in which the four-cycle volumes are stabilized nonperturbatively, the 
b axion receives a mass in a stabilized vacuum. 

At this stage the mass-squared ml of b is proportional to the vacuum energy and hence is 
negative in the supersymmetric AdS minimum. The minimum of the potential will be the final 
point of the inflationary dynamics, and hence we would like it to have a very small positive 
cosmological constant to be consistent with the current accelerated expansion of the universe. 
Thus, we need to include an uplifting term. In the uplifted minimum, ml oc Vds > 0. This 
relation is the origin of the eta problem that was found in [S] choosing b as an axion: for a 
generic uphftingj 1 ^] V" (b) ~ V(b), so that rj ~ Oil) and slow roll inflation does not take place. 
This is completely analogous to the eta problem of D-brane inflation found in [39J and can 
be intuitively understood in the same way. Here we will take b = as the stabilized value _ 
of b and concentrate on c as a candidate inflaton. 

Let us now turn to consider c, which does not appear in the Kahler potential or superpo- 
tential at any order in perturbation theory. To assess c as an inflaton, one should determine 
the leading nonperturbative effects, either in the superpotential or in the Kahler potential, 
that do introduce a potential for c, i.e. one should identify the leading breaking of the shift 
symmetry. Euclidean D3-branes carrying vanishing Dl-brane charge do not induce a poten- 
tial for c, but Euclidean D3-branes supporting worldvolume fluxes (and hence nonvanishing 
Dl-brane charge) give rise to a dependence on c, via the Chern-Simons coupling J F 2 A C%. 
As observed in [5], it follows that when the Kahler moduli are stabilized by Euclidean D3- 
branes, c receives a mass in the stabilized vacuum: one must sum over Euclidean D-brane 
contributions to the superpotential, including summing over the amount n — J F 2 of magne- 
tization, and this generically introduces an eta problem for c. The solution, as explained in 
[5], is to stabilize the Kahler moduli via gaugino condensation on D7-branes, which leads to 
an exponentially smaller (and hence negligible) mass for c. 



14 Notice that the proportionality constant in oc Vds depends on the volume-dependence of the uplifting 
term and could be made small for particular choices of the latter as proposed in [35] . 

15 It is easy to check that b — is still the stabilized value after the inclusion of nonperturbative corrections 



to the Kahler potential, cf. {6.5 
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5.3 Axion decay constants in string theory 

We now turn to the important task of expressing the axion decay constant, /, in terms 



of the data of a compactification. As we reviewed in §5.1[ the decay constant of an axion 
Ci = c(x)lj/(2tt) is given by 



P 



M 2 
p 



J uj A *u> 
A8n 2 V E [ (2vr) 6 a /3 



9s 



(5.18) 



so that the primary task is to compute the norm J uj A *uj. (This problem has been studied 
in a wide range of examples in [40J . ) We will first recall, in £5.3.1, how to express the axion 
kinetic term, and hence also the axion decay constant, in terms of j\f = 1 data. This will 
lead us to a simple expression for the decay constant in terms of intersection numbers of the 
Calabi-Yau. We will then propose a class of models in which the decay constant is rather 
small, motivated by the fact that with other parameters held fixed, decreasing / increases 



the amplitude of the resonant non-Gaussianity. Next, in £5.3.2, we will present a concrete 
example that illustrates the geometry of a configuration that leads to small /. 



5.3.1 Decay constants in terms of j\f = 1 data 



In £5.2.1 we have reviewed, following [35], the four-dimensional j\f = 1 description of Type IIB 



03-07 orientifolds. The multiplets relevant for us are the orient if old-odd chiral multiplets 
G a and the orientifold-even chiral multiplets T a . The tree- level Kahler potential given in 



(5.9) determines the kinetic terms for G a and hence the decay constants of the axions b a and 
c a . First let us notice that the Kahler metric in the space of the chiral multiplets T a and 
G a factorizes in two blocks, K^ a f and K G a G b. The reason is that off-diagonal terms such as 
K Ta Qa are proportional to intersection numbers c a p a with one odd index and two even indices, 
which are forbidden by the orientifold action [22]. We are interested in one particular mode 
from among the G a , which we will denote by £_ is then the orientifold-odd two-cycle 
that supports our candidate inflaton c~. We now choose a basis for G a such that K G a G b is 
block diagonal with a 1 x 1 block K G - G ~. The kinetic term for c~ is then given by 



M 2 p K G -a- 



(2vr) : 



{do') 2 cM*K G - G -\dG- 



where 



K. 



G~G~ 



d 2 K(G, T) 
dG-dG- 



4V F 



(5.19) 



(5.20) 



and we used 



2t q + g s c abc Im G b Im G c 



(5.21) 
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Hence we can express the decay constant of the axion c as 

P 9s c a ^v a 



Ml 8tt 2 V e 

As promised, we have expressed the norm J uj A *uj in terms of the intersection numbers 

J uj A *uj 2 



(5.22) 



(27r)V 3 3 



c a __v a . (5.23) 



In £6.4 we will discuss the constraints that follow from the result (5.22). First, in the following 



subsection we provide some geometrical intuition for (5.22) 



5.3.2 An example: a complex plane of fixed points 

An instructive example arises from considering an orbifold that is locally C 2 /Z 2 x C, i.e. an 
Eguchi-Hanson space fibered over a base S of complex dimension one. Let uj be the two-form 
dual to the blowup cycle of the orbifold, and let E be the two-manifold of fixed points, i. e. the 
base over which the Eguchi-Hanson space is fiberedj^] We are interested in the decay constant 
of the axion C 2 = ^-c(x)uj, so we must compute J uj A *uj. In the local approximation, this is 
straightforward, as we shall see. However, far from the fixed-point locus, the fiber may deviate 
substantially from the Eguchi-Hanson geometry, in a complicated and model-dependent way, 
and moreover the fixed-point locus E may be embedded in the compact space in a nontrivial 
manner. Happily, the integral f uj A *uj has its primary support near the fixed-point locus, 
where the local approximation is excellent. 

We recall, following the useful summary in Appendix B of jH], that the Eguchi-Hanson 
space has a unique homology two-cycle of radius a/2, where r = a defines the location of the 
coordinate singularity; here r is the standard radial coordinate. The two-form uj corresponding 
to this cycle may be written 

q 2 / dv dv I \ 

uj = 2-nol— — A dip + cos 9— A d<p + -sin Q dQ A d<f)) (5.24) 
r 1 \ r r 2 / 

in terms of r and the angular coordinates ip,8,<j). By observing that *^uj = —uj and that 
/ uj A uj = -(2tt)V 2 /2, one finds 

J EH UJA*4UJ _ 1 

(2n)W ~ 2 ' 

Clearly, given the form of uj, this integral has its support in a region a < r < few x a. 
This justifies the local approximation as long as the compact space has a radius that is large 



16 Concretely, we are imagining that £ extends into a warped throat region, and that an NS5-brane wraps 
the blowup cycle at a particular location in the throat. The warping is invoked in order to suppress the energy 
density of the wrapped NS5-brane. See [5] for further details, and for an example of a suitable orbifold action 
in a Klebanov-Strassler throat. 
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compared to a. Next, we observe that 



j^ A *6^ = f EH d*xuAnu J sv ^ = 1 , . 
(2vr) 6 a' 3 (27r) 4 a /2 (2tt)V 2 1 } 



(5.26) 



Substituting this in (5.18), we recover the parametric scaling of £5.3.1 



6 Microscopic Constraints 

We now turn to determining the ranges of our phenomenological parameters that are allowed 
in a consistent and computable microphysical model. 

Let us first remark that, as usual in string theory model building, computability imposes 
stringent constraints on the compactification parameters. Because large-field inflation in- 
volves substantial energy densities and requires correspondingly steep moduli barriers, the 
compact space needs to be reasonably small, so that the Kaluza-Klein scale and the (neces- 
sarily lower) scale of moduli masses can be large enough to prevent runaway moduli evolution. 
Clearly, one must then carefully check that the compactification is still large enough for the 
supergravity approximation to be valid; furthermore, backreaction of the inflationary energy 
on the compact space is a serious issue, particularly when this space is not large in string 
units. Incorporating these requirements then leads to severe restrictions on the allowed values 
of the decay constant /. 

We will begin in §6.1| by considering the constraints from computability, then give, in §6.2[ 
a qualitative description of the constraints from backreaction, deferring details to appendix B. 



Next, in §6.3[ we will verify that a two-derivative action suffices to describe this system. This 
is not obvious, as rapid oscillations in the potential could enhance the importance of generic 
higher-derivative terms; however, we will show that the specific terms emerging from string 



theory are negligible in our solution. We then apply these constraints in £6.4 to determine 



the range of the decay constant /. Finally, we estimate the size bf of the modulations; as 



this is rather model-dependent, in §6.5[ we will restrict our attention to a specific example in 
which a periodic contribution is generated by Euclidean Dl-brane corrections to the Kahler 
potential. 



6.1 Constraints from computability 

In this subsection we will list several constraints coming from the consistency of the string 
theory setup. We will first require the validity of the string and a' perturbation expansions, 
and the validity of neglecting higher-order corrections to the nonperturbative superpotential, 
and then we will require that the inflaton potential does not destabilize the compactification. 

First of all we require the validity of string perturbation theory, i. e. we require j s < 1. We 
must also ensure the validity of the a' expansion. To do this including a reasonable estimate 
of numerical factors such as 27r, it is convenient to use worldsheet instantons as a proxy for 
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perturbative a' corrections, because the normalization is easily determined. To get the correct 
coefficient, we start from the string-frame ten-dimensional metric g string and impose that the 
worldsheet instanton action obeys e~ Sws < e~ 2 , or 



l ~, J V '9 string > 2 , (6.1) 



2ira 

which using g string = gEinstein\fg~s is converted to Einstein frame 

2<&fj = ^ =, (6.2, 



and we used that L a cj' 3 = (27r) 2 a / (5 Q , /3 . 

As we invoked nonperturbative corrections to the superpotential, we must also require 
that any further superpotential corrections, e.g. from multi-instantons, are negligible. For 
this purpose it suffices to impose 

e -a a T a < e - 2 ^ 1 ^ r Q > — . (6.3) 

7T 

Additional constraints come from the moduli stabilization process. To use the single-field 
inflationary analysis we have developed in §|2] and £|3j we need to require that the uplifted 
minimum is only slightly perturbed by the inflationary dynamics. In particular, the linear 
potential that we have represented as fj?(f> actually depends on the compactification volume, 
and hence shifts the minimized value of the volume. In four-dimensional Einstein frame, the 
leading term in the inflaton potential is 



(MY 



^,Vi)« [y^J M (6.4) 

where (Ve) is the expectation value of the volume. To ensure that the resulting contribution 
to the potential for the volume is unimportant, we will insist that the inflaton potential 
induced by the NS5-brane, V(<p), is smaller than the moduli potential U mo d- 
At the supersymmetric minimum we have 



V A dS = ~^^2~ ■ ( 6 ' 5 ) 



Without specifying the details of the uplifting mechanism, we assume that an uplifting to 
a small and positive cosmological constant is possible, and that the height of the potential 
barrier U mo d that separates the uplifted minimum from decompactification is of the same 
order as U mo d ~ |Va<2s|. Now, the COBE normalization tells us that 

V(<Pcmb) = e ( 0.027M P ) 4 ~ 2.4 • !Q- 9 M p 4 . (6.6) 



34 



Hence we obtain the constraint 

9s Z\W\ 2 
2 VI 



w, 



AdS\ 



~ W mod > 2.4 • 10" y M o 4 . 



(6.7) 



To extract a useful form of the above constraints, let us substitute for W the solution of 



any of the equations (5.11) 



W = +\A a \a a e 



2V F 



(6.8) 



with no sum over a. We will also assume \A a \ ~ 1 (see [5] for a discussion of this point). 
After some manipulations we find 



r Q <-^log( A,, I0 



(6.9) 



again with no summation over a. Finally, we should limit the number of D7-branes in each 
stack; although there plausibly exist examples with N a quite large, we will impose N a < 50. 
This gives us 



r Q < 73 - 8 log 



2g s 



(6.10) 



We notice that v a {ix^/g~ s ) > 1 was the condition in (6.2) that enabled us to neglect a' cor- 
rections, so that as long as g s ^ 0.5 the second term on the right hand side of (6.10) is 
negative. 



6.2 Constraints from backreaction on the geometry 

Another important constraint comes from the requirement that the backreaction of the infla- 
tionary energy density on the compact space is small. In this section we will give a qualitative 
description of the problem and will briefly sketch a model-building solution; the interested 
reader is referred to appendix B for a more complete treatment. 

At the time that the CMB perturbations are produced, the inflaton has a large vev in 
Planck units, <fi ~ 11M P , corresponding to a configuration of the two-form potential threading 
the two-cycle £_ of the form 

/ ^ = iV. = ^ 7 »l (6.11) 

(27T)W ,/ s _ 27T/ 

In the absence of an NS5-brane wrapping E, there would be no energy stored in this configu- 
ration, as C2 enjoys a shift symmetry. However, inflation is driven by the substantial energy 
stored in this system by the Born-Infeld action of the wrapped NS5-brane. Moreover, there is 
a corresponding D3-brane charge induced by the Chern-Simons coupling J C2AC4. Note that 
the net induced D3-brane charge in the total compactification is zero, as required by Gauss's 
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law, because we have arranged for an additional, tadpole-canceling NS5-brane that wraps a 
distant cycle homologous to £_, but does so with opposite orientation. Therefore, the 
Chern-Simons coupling induces a dipole configuration of D3-brane charge, with F§ flux lines 
stretching from XL to £'_. 

It is essential to ensure that the inflationary energy, which is effectively localized in the 
compact space in the vicinity of the wrapped NS5-brane, does not substantially correct the 
remainder of the compact geometry. Heuristically, one can imagine that the increased tension 
of the NS5-brane, as well as the induced charge, is represented by N w D3-branes dissolved in 
the NS5-brane. We must therefore estimate the effect of N w D3-branes in a warped throat 
(recall that we have situated each wrapped NS5-brane in a warped region in order to suppress 
its energy density below the string scale, as required e.g. by the COBE normalization). 
Clearly, this backreaction will be reasonably small if N w <C N, with iV the D3-brane charge 
of the background throat. 

However, we must be careful about the effect of even a modest distortion of the geometry 
on the moduli stabilization and therefore on the four-dimensional potential. Let us first recall 
that in scenarios of D3-brane inflation in nonperturbatively-stabilized vacua, even a single D3- 
brane moving slowly in a throat can affect the warp factor, and correspondingly the warped 
volumes of four-cycles bearing nonperturbative effects, to such a degree that this interaction 
is the leading contribution to the inflaton potential [121 H3] . 

This sensitivity originates in two facts: first, D3-branes perturb the warped metric in a 
manner that is not suppressed by the background warp factor at the location of the D3- 
branes, because D3-branes are BPS with respect to a throat generated by D3-brane charge, 
and hence their contributions to the metric may simply be superposed on the background. 
Second, nonperturbative effects on a four-cycle are exponentially sensitive to changes in the 
four-cycle volume. Both these facts appear threatening for a situation such as ours in which 
the moduli are stabilized nonperturbatively and substantial D3-brane charge is induced in a 
throat: one can anticipate that as inflation proceeds and the D3-brane charge diminishes, the 
four-cycle volume changes, leading to an unanticipated, and possibly steep, contribution to 
the inflaton potential. 

To understand this concretely, we will first consider a simpler system: an anti-D3-brane 
in a warped throat generated by N D3-branes, or equivalently a warped throat generated by 
N — 1 D3-branes, together with a brane-antibrane pair. Furthermore, from the result of [S] 
one learns that at long distances, the effect of the brane-antibrane pair on the supergravity 
solution is strongly suppressed by the warp factor at the location of the pair, i.e. at the tip of 
the throat. In contrast, the effects of D3-branes are not suppressed in this manner. Therefore, 
for the purpose of computing perturbations to the bulk compact space, we may replace an 
anti-D3-brane in a warped throat generated by iV D3-branes with a warped throat generated 
by iV — 1 D3-branes, up to exponentially small corrections. 

Equipped with this approximation, we may represent the configuration of interest as fol- 
lows: two warped throats, carrying the charge of Ni, N 2 D3-branes respectively, are perturbed 
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to Ni + N w , N 2 — N w by the inclusion of the NS5-brane in (say) the first throat, and the anti- 
NS5-brane in the second throat. Here we are ignoring the warping-suppressed correction 
indicated above, and we are approximating the NS5-branes by the D3-brane charge and ten- 
sion that they carry, which is an excellent approximation for N w ^> 1. Other effects due to the 
NS5-brane that do not depend on its induced D3-brane charge, i.e. on its world- volume flux, 
are independent of the inflaton and hence do not correct its potential. One can now easily see 
that the volume of a four-cycle at a generic location in the compact space will be corrected 
by the inclusion of the NS5-branes. If the four-cycle happens to enter one or both throats, 
the change in the volume is easily computed, and is seen to be substantial (cf. appendix B). 

To control this problem, we situate the NS5-brane and the anti-NS5-brane, together with 
the family of homologous cycles connecting them, in a single warped region. The idea is 
that from the bulk of the compact space, the NS5-brane configuration will appear to be a 
distant dipole whose net effect, integrated over a four-cycle, averages out to be small. This 
setup allows us to parametrically suppress the backreaction by a small factor given by the 
ratio of the dipole length, i. e. the distance between two NS5-branes, to the distance between 
the NS5-branes and the four-cycle in question. This small factor comes in addition to the 
suppression by the small ratio N w /N^\ 

In appendix [B] we give more details about the above setup. We show, through two explicit 
models of increasing complexity, the robustness of the above suppression mechanism. 



6.3 Constraints from higher-derivative terms 

The analysis presented thus far has used the two- derivative action, which is an approximation 
with a limited range of validity. In general, one expects an infinite series of higher- derivative 
terms, possibly including multiple derivatives as well as powers of the first derivative. Our 
background solution involves rapid oscillations, so it is reasonable to ask whether these high 
frequencies enhance the role of higher- derivative terms and render the two-derivative approxi- 
mation invalid. To check this, one should evaluate the higher-derivative terms on the solution 
and compare to the two- derivative action. We will now show that the two-derivative ap- 
proximation is valid in the scalar sector; analogous considerations apply to the gravitational 
action. 

Rather than write down the most general higher-derivative corrections to the scalar sector, 
we give here the terms that end up being present in the string theory examples. In string the- 
ory, we can directly compute the leading higher-derivative terms in the action for b, extending 
the result to c using S-duality. To get the leading terms, one considers the a' 3 corrections to 
the effective action due to Gross and Sloan jl5] (at the four-point level) and Kehagias and 
Partouche [IB] (up to the eight-point level). These corrections are of the same lineage as the 

17 A further suppression can be achieved with a carefully-chosen embedding of the four-cycle, e.g. one that 
is symmetric with respect to the two NS5-branes. However, this requires fine-tuning, whereas the dipole 
suppression on which we have focused is parametric. 
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famous Riemann 4 term, but involve NS-NS three-form flux. This yields corrections to the 
axion kinetic terms. Following |46j, the ten-dimensional Einstein-frame action including the 
leading (a' 3 ) corrections is 



SlOD,E 
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1 



(2yr) 7 a 



d w Xy/gE I R E - 



1 



I2g. 



it ttKMN 



+ 



R 



PQ 



MN 



R 



PQ 



C(3) 
3 x 2 6 



+ 2^ W [M^N] ~ ids H [M h N]C ~ • ■ • 



//, :i 2 n ,: '7." • .(6.12) 



Q] 



MN 



(6.13) 



and the square brackets are defined without the combinatorial factor 1/2 in front. Hence, 
the terms that are relevant for our axion at order a' 3 are proportional to if 8 and (VH) A . To 
estimate the importance of these terms, we will consider a special case in which the internal 
space is a T 2 x T 4 , with the NS-NS two-form field only along the T 2 directions 8 and 9, i.e. 
B 89 = —B 9S = b. Furthermore, since the background dynamics involves large frequencies 
but not large spatial gradients, we are primarily interested in terms containing only time 
derivatives, and can therefore take b to be homogeneous in the noncompact spatial directions. 
In this special case, making use of (2.13) in [15], and using S-duality to determine the action 
for c from that for 6, we find that after dimensional reduction the corrected action for c is 
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cf to make the kinetic term canonical, yielding the action in terms of 
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where we can now calculate the scale of the higher derivative terms Mj and Mjj to be 

i/ 



Mi = M, 



9s 



C(3) 



(6.16) 



and 



M: 



M, 



3/16 

9* 



3 9 10\ V8 



TT 2 



c(3); 



(6.17) 



To determine whether these higher-derivative terms will become important, we compute the 
dimensionless quantity m jj i , where u — j (3.5) is the physical frequency of oscillations; we 
obtain 



38 



Mu " V IO- 3 J V0.2/ V 120 ' 



For the ranges of / and V# that will be of interest to us (cf. the higher-derivative terms 
are not important and our two-derivative approximation is justified. 

6.4 Constraints on the axion decay constant 

In this section, we discuss direct constraints on the axion decay constant /. We first recall a 
rather general (conjectured) upper bound / < M p [47J, and we then describe and incorporate 
a novel lower bound, specific to our setup, that arises from combining the requirements that 
a' perturbation theory should be valid and that the inflationary energy should not drive 
decompactification. 

Despite many attempts, at the time of writing there is no known, controllable string theory 
construction that provides / > M p . In particular, the authors of jH] have scanned several 
classes of string theory models and found sub-Planckian axion decay constants in every case. 
However, this upper bound on / is of relatively little importance for the phenomenological 
signatures we are considering in this paper. 

On the other hand, a potential lower bound on / is of considerable importance for our 
analysis. Considering oscillations in the CMB spectrum, in the regime / -C M p one can 
easily find models that range from being observationally excluded to giving undetectably 
small modifications, depending on the amplitude of the ripples in the inflationary potential. 
Furthermore, the resonant non-Gaussianity becomes large only for small / (e.g. we will find 
that / < 3 • 10~ 3 is a necessary condition to give a reasonable prospect of detectability). 
Hence we will move on to consider possible lower bounds on /. 

As discussed in [5] and in the preceding section, a direct lower bound on / comes from 
the requirement of small backreaction. In particular, the radius of curvature induced by 
the energy localized on the wrapped NS5-brane should be smaller than the smallest radius 
of curvature R± in a direction transverse to the NS5-brane in the compactification. This 
requires 

Rjx f 

4ng s M p " R\X 

where we have defined X = Vol(X5)/7r 3 , with X§ the base of the cone forming the warped 
throat. We remark that X < 1, as S 5 is the Sasaki-Einstein manifold with the largest volume, 
in the sense defined above. We can estimate R± as being comparable to the AdS radius R of 
the throat containing the NS5-brane. Given that the volume^ V of the Calabi-Yau has to be 



N ™^l±r (6-20) 
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We always refer to the warped volume, calculated with the whole warped metric. 
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larger than the volume of any throat it includes, one finds that 
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where for simplicity we have assumed that the UV cutoff of the throat is at r ~ R where the 
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warp factor becomes of order unity. Putting together (|6.20|) and (|6.21|), we find 

/ 



Mr, 



> 



0.09 



^^9s 

XV3y2/3 XV3V2/3 xVsy 273 



(6.22) 



Although the above constraint substantially restricts our parameter space, an even stronger 



constraint comes from demanding the validity of a' perturbation theory: using (5.22) for / 



and combining this with the lower bound on two-cycle volumes given in (6.2), we obtain 
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where we have assumed that c a > 1. (6.23) turns out to give the strongest microphysical 



lower bound on /. An upper bound is harder to determine from this formula. Assuming 

) p 



again that c a > 1, assuming that no precise cancellations occur, and using (6.2), we find 
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(6.24) 



6.5 Constraints on the amplitude of the modulations 



So far we have seen that with the Kahler potential and superpotential given in (5.9) and 



(5.10), the axion c persists as a flat direction after moduli stabilization, 19 As explained in [5] 



the presence of an NS5-brane wrapping the two-cycle that defines c introduces a monodromy 



and results, for large c, in the linear potential in (5.16). In this section we will consider 



further nonperturbative corrections that will in general induce small modulations of this 
linear potential. These are precisely the modulations whose phenomenology we have studied 
in the first part of this paper. 

Nonperturbative corrections could appear both in the Kahler potential and in the superpo- 
tential. We focus on the first possibility and comment at the end of this section on the second. 
Consider the type IIB orientifolds with 03-planes and 07-planes. As we have remarked, the 
RR two-form C2 is odd under the orientifold projection and therefore a four-dimensional ax- 
ion that survives projection comes from integrating C2 over an odd two-cycle v~ . Such an 



odd cycle can be thought of as v~ 



', where v 1 and v 2 represent two two-cycles in 



the parent Calabi-Yau manifold that are mapped into each other by the orientifold action. 



19 As we have remarked, the axion b has its flat direction lifted by nonperturbative stabilization of the 
Kahler moduli. 
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Now consider a Euclidean Dl-brane wrapping the even cycle v + = v 1 +v 2 . Such an instanton 
feels the local M = 1 supersymmetry of the orientifolded theory, and it breaks this super- 
symmetry completely]^] Hence this is a non-BPS instanton with four universal fermionic zero 
modes, namely the goldstini of the broken M = 1 supersymmetry. If the Euclidean Dl-brane 
wraps a minimum-volume cycle in the homology class v + then it has the right total number 
of fermionic zero modes (four) to contribute to a D-term and in particular to the Kahler 
potential. 

More specifically, in [49J it was argued that nonperturbative contributions from worldsheet 
instantons and their SX(2,Z) images, Euclidean (p,q) strings, give rise to corrections to the 
prepotential of the M = 2 theory of the parent Calabi-Yau compactification. Such corrections 
are most naturally expressed inside the logarithm of the Kahler potential, 

K = -2\og[V E + g(G,G)} , (6.25) 

where g is an appropriate function. Invariance under SL(2, Z), or more generally under a 
subgroup T C SL(2, Z), is naturally achieved if g is the sum of some individual correction g 
over an orbit of T. 

At the time of writing, the nonperturbative correction g is not known explicitly, but 
a modular-invariant result has been conjectured in [50]. Inspired by the structure of this 
result (which we will not reproduce here), we will make a simple educated guess based on 
the following criteria: the non-perturbative correction should go to zero exponentially for 
large two-cycle volume v + ; it should break the continuous shift-symmetry of c to a discrete 
shift-symmetry c — > c + 2ir; and it should be invariant under whatever discrete subgroup 
T C SL(2, Z) of the ten-dimensional SL(2, Z) symmetry is preserved by the compactification. 
The subgroup T may well be trivial, and we will assume this for simplicity; note, however, 
that one can plausibly obtain a more constrained result when some or all of the symmetry 
is preserved, as in [50J. Moreover, notice that along the orbits of T, the instanton action 
generally increases compared to that of a single worldsheet instanton or Euclidean Dl-brane; 
thus, when the volume v + is not too small, only a few terms make an important contribution, 
with the remainder enjoying further exponential suppression. 

A reasonable guess satisfying these criteria, for T trivial, is 

K = —2 log [Ve + e~ SED1 cos(c)] = -2 log 

In light of this corrected Kahler potential, we should revisit the moduli stabilization before 
proceeding to calculate the size b of the periodic contribution to the scalar potential. 

20 To see this, note that (cf. [38] ) the instanton action depends on a two-cycle volume, but the proper Kahler 
coordinates arc four-cycle volumes. Therefore, the instanton action cannot be holomorphic, so the instanton 
cannot contribute to a superpotential, and must instead be non-BPS. 



Ve + e ^ cos(c) 



(6.26) 
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We begin by noticing the following implication 
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which can be verified by direct computation. This allows us to use the F-flatness condition 
to find the minimum in the T Q -directions even when one of the F-terms, namely D G _W, does 
not vanish. Equipped with this knowledge we repeat, mutatis mutandis, the steps of §(5l 



First, the phases of the T a are stabilized as in (5.13), with k being odd as explained below 



(5.14). The reason is that the sign of dr a K is not changed by the small nonperturbative 



correction e Sedi . Second, ImG is again stabilized at 0. Given (6.27), the equation one needs 
to solve is D G aW = 0, which reduces to 
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where we made use of (5.6) to perform the derivative on c. Since Ve and v + only depend on 



ImG implicitly as in (5.21), we can take the imaginary part of (6.28), 
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But from (5.21) we know that 

c a /3 1 (di mG ^v /3 )v' Y = g s c aa -lmG a 



(6.30) 



which means that ImG" = (for every a) is a solution to (6.29). The real part of (6.28) 



is nonvanishing and of order e Sedi . Again because of (6.27), the minimization in the r Q 
is obtained by imposing D Ta W = 0. These equations depend on the inflaton c, appearing 



explicitly in (6.26), and hence the minimum in the T a directions will be a function of c. 



Integrating out the T a leads to a contribution in the effective potential V[T(c), c] for c which 
is of the same order as the contribution coming from the explicit c-dependence in the Kahler 
potential. Therefore this effect cannot be neglected. To take it into account, we solve the 
Dr a W = equations perturbatively in e~ SED1 . 

We define the coefficients of the minimum in the r a directions in a perturbative expansion 
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and so on for all other variables. The zeroth-order equations are 



(D Ta W) {0) = (d Ta W) {0) + W (0) (d Ta K) m = 



(6.32) 
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which can be solved numerically once the model is specified. The first-order equations are 



(D Ta W) {1) = (d Ta W) {1) + W (1) (d Ta K) {0) + W {0) {d Ta K) {1) = 



(6.33) 



which again can be solved numerically using the solutions of (|6.32|). We turn now to estimate 

W (l) 



the parameter b defined in (2.1). One finds 
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where we have defined U m od as the moduli stabilization barrier at zeroth order in e~ 
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More explicitly, using ( |6.26[ ) and ( |5.10[ ), 
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where the first line is valid for any a and the second line (obtained using (6.33)) shows that 
the expression for b is independent of a. Notice that dT a v~ 

Capjidr^v 13 



d To r a 



\d Ta v + is given implicitly by 

(6.37) 

Some comments are in order. The size of the ripples in the potential is proportional to the 
ratio of the moduli stabilization barrier to the scale of inflation, which has to be large for the 
self-consistency of the estimate. We have used the value of the potential at the would-be AdS 
minimum to estimate the moduli stabilization barrier once an uplifting term is included. Due 
to the exponential suppression e~ SED1 , the size of bf is extremely sensitive to g s and v + . 

An upper bound can be derived from (6.36) using the following considerations. In the 
KKLT construction, perturbative corrections to the Kahler potential can be neglected as long 
as Wo 1, and generically W ~ Wq. For larger values of Wq, perturbative corrections have 
to be included, as in the large volume scenario [51J. In the present work, we focused on 
the former setup and we leave an investigation of the latter for the future. The exponential 
suppression in (6.36) can be bounded by (6.2). Finally, we denote the model-dependent term 
in square brackets in (6.36) by c . Putting things together leads to the bound 
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Even imposing all the model-independent constraints we have described in the previous 
sections, one can still have bf > 10~ 4 , which, as shown in ^4j is roughly the upper bound 
imposed by measurements of the scalar power spectrum. Therefore, in certain parameter 
ranges the primary constraint on modulations of the potential comes from the data, not from 
microphysics. 
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7 Combined Theoretical and Observational Constraints 



We now summarize our results, combining the observational constraints from ^4] with the the- 
oretical constraints from §[6j As an aid to the reader, we will now briefly recall the qualitative 
properties of those results. 

Axion monodromy can produce characteristic signatures in the CMB: the oscillations in 
the axion potential generated by nonperturbative effects source resonant contributions to the 
scalar power spectrum and bispectrum. The amplitude and frequency of the oscillations in 
the potential can therefore be bounded by comparison to observations. We recall from §[4] 
that the observational constraints take the form of exclusion contours in the space of the 
phenomenological parameters, after marginalization over additional model parameters that 
have important degeneracies with those displayed. For convenience, we have chosen to display 



constraints in terms of the parameters / and bf defined in (2.1), marginalizing over the phase 
A(f) and over f^/i 2 . 

The first new step is to combine the exclusion contours based on the temperature two- 
point function with estimates of constraints from the three-point function. Based on the 



rough estimates described in §3.4[ we present, in figure [6j three contours at f res = 200,20,2, 
with the expectation that the gray region (f res > 200) might plausibly be excluded, while 
the colored, lighter regions (20 < f res < 200 and 2 < f res < 20, respectively) are possibly 
within detectability. A careful study of the constraints on resonant non-Gaussianity would 
be a worthwhile topic for future research. 

Next, we recall that in §[6j we found that the requirements of consistency and computability 
in the string compactifications giving rise to axion monodromy models led to constraints on the 
parameters / and bf. Let us remark that as these constraints are not rooted in deep principles 
of string theory or of quantum field theory, but rather originate in practical limitations in our 
present ability to construct computable models, they may well be loosened in further work. As 
such, the theoretical constraints we present here should be understood as designating included 
rather than excluded regions: in contrast to experimental contours, theoretical contours of 
this sort may expand rather than contract given improved understanding. 

Because the parameter b measures the amplitude of a nonperturbative effect, it is expo- 
nential in the natural input parameters, and can therefore be made small without substantial 
fine-tuning. We therefore do not present a lower bound on b. However, we found the theoret- 
ical upper bound ( |6.38f 



21 



bf < 2c • 10 9 M„ 4 ^-e- 2/9s , (7.1) 



V 



with a model-dependent constant Co that can be estimated in explicit examples, and which 
we find to be typically of order 10~ 2 . 



21 We stress that this 'bound' is not universal and depends on the assumptions enumerated in ^6] We include 
it here as a representative example of the constraints that arise in particular scenarios. 
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Figure 6: We show the (one- and two-sigma) likelihood contours for the temperature two- 
point function together with three contours that characterize the amplitude of the three point 
function, for f res = 200, 20, 2. 



Next, we obtained a lower bound for / in (6.23) A precise upper bound, however 



is highly model-dependent. We estimate an upper limit by assuming that the intersection 



numbers are of order one^jand that no precise cancellations occur. From (6.23) and (6.24) 
the complement of the theoretically excluded range for / is then 
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(2tt) 3 / 2 , 

Notice that the theoretical constraints depend mainly on two quantities: the string cou- 
pling g s and the volume Ve of the compactification. The former appears in the exponential 
suppression of the nonperturbative effect generating the modulations of the potential. Hence, 
g s < 0.1 suppresses any possible signature of the modulations. For g s > 0.1, there is always 
a theoretically allowed region in which the oscillations in the inflaton potential lead to ob- 
servable ripples in the two-point function of the CMB. On the other hand, the size of the 
non-Gaussianity depends critically on Ve as well. Assuming g s > 0.1, larger Ve allows for a 



larger range of / and therefore larger non-Gaussianity (see (7.2)). A way to quantify this is 



to use the estimate obtained in £ 3.4 
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is weaker than (6.23) 



22 Thc constraint from the backreaction described in {t 
23 Larger intersection numbers are an interesting possibility that we will not investigate here. 
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Figure 7: We superimpose the theoretical constraints, summarized in (7.1) and (7.2), on 
the constraints imposed by observations, which are shown in figure [6j The orange overlay 
indicates regions of the parameter space that are difficult to reach in the class of models 
considered in the present work. The theoretical constraints are shown for g s = 0.12,0.5 and 
V E = 100, 900. 
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Figure 8: The blue dot represents the explicit numerical example presented in appendix [DJ It 
represents a case in which upcoming experiments could detect the signatures of modulations 
in both the two-point function and the three-point function. 



and the lower bound in (7.2). The result is 



4 /5 /m-4\ 4/5 



bf 

We now combine the theoretical and observational constraints, presenting them in the 
plane {log(/),log(6)}. We choose as boundaries 10~ 4 < f <^ M p and 10~ 4 < b <C 1 based on 
the following considerations. The number of oscillations per e-folding is roughly 10~ 2 M p / f . 
Hence for / > 0.1M p there is less than one oscillation in the whole range of scales probed by 
the CMB, and the signal from modulations becomes degenerate with the overall amplitude. 
Furthermore, in £j3j we systematically used the expansion b <C 1, where 6 = 1 divides mono- 
tonic from non-monotonic potentials. Finally, the lower boundaries 10~ 4 < / and 10~ 4 < b 
exclude regions that are relatively uninteresting in the present context: smaller values of b 
lead to an unobservably small signal, while smaller values of / are rather difficult to obtain 
in the class of string theory constructions we considered. In the {log(/), log(6)} plane, the 
theoretically allowed region looks like an interval in /, whose size is determined by Ve, with 



an upper cut effectively determined by g s as in (7.1). 

Finally, in figure[8]we show where a particular numerical toy example, with specific choices 
of the intersection numbers, lies in the {log(/), log(6)} plane. 
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8 Conclusions 



The goal of this investigation was to characterize the predictions of axion monodromy inflation 
for the CMB temperature anisotropies. Nonperturbative effects in these models generically 
introduce sinusoidal modulations of the inflaton potential, which in turn lead to resonantly- 
enhanced modulations of the scalar spectrum and bispectrum. 

We have provided a simple analytic result for the modulated scalar power spectrum in this 
class of models. We also presented an alternative derivation in terms of episodes of particle 
production driven by resonance between a mode inside the horizon and the driving force of 
the oscillatory background evolution. We then determined in detail the constraints that the 
five-year WMAP data places on models with modulations of this sort. 

Next, after reviewing the realization of axion monodromy inflation in string theory, we 
performed a comprehensive study of the parameter constraints implied by the requirements 
of microphysical consistency and computability. The resulting allowed parameter regions are 
very plausibly realizable in sensible string theory constructions. 

We also identified a new contribution to the inflaton potential in axion monodromy infla- 
tion: the backreaction of the inflationary energy on the compact space can source an important 
correction to the potential by correcting the volumes of four-cycles and hence affecting the 
scale of nonperturbative moduli-stabilizing effects. We then presented a model-building solu- 
tion to this problem, in which the NS5-brane and anti-NS5-brane driving inflation are in the 
same warped region, or more generally are distant from the four-cycles of interest. 

Finally, we combined the observational and theoretical constraints, in order to ascertain 
whether detectable modulations of the scalar spectrum and/or bispectrum are possible, con- 
sistent with current observational bounds and known theoretical restrictions. Our conclusion 
is that both sorts of modulations are possible, and in fact in many cases the strongest bound 
on the amplitude of the modulations comes from data, not from microphysics. Moreover, 
even though observational limits on the amplitude and frequency of modulations in the scalar 
power spectrum provide a strong constraint on the parameter space of axion monodromy 
models, and even though microphysical constraints sharply restrict the allowed frequencies, 
detectably-large non-Gaussianity can indeed be produced in a class of controllable models. 
Such models enjoy three nontrivial signatures: detectable tensors with r 0.07, a modulated 
scalar power spectrum, and resonant non-Gaussianity. 

Let us remark that even in the absence of non-Gaussianity, this class of models is eminently 
testable: axion monodromy inflation unambiguously predicts a large tensor signal, and the 
parameters of the models are already strongly constrained by limits on modulations in the 
scalar power spectrum. 

There are several interesting directions for future work. First, we have not analyzed 
the constraints on the model from the three-point function; more generally, understanding 
the prospects for constraining or detecting resonant non-Gaussianity is an important task. 
Moreover, it would be instructive to construct an explicit model in which the many theoretical 
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constraints we have checked can be combined in a coordinated way. In addition, it would be 
interesting to determine whether chain inflation can be realized in this context. 

It is intriguing that the modulated power spectrum we have found is very similar in form 
to that proposed in the context of modifications of the initial state, as in e.g. [22 | 1231 l8| 191 ITU] . 
In light of our calculation of the power spectrum in terms of particle production (£ |3.3 ), this 
coincidence is not entirely surprising: the driving force of the oscillating background eventually 
generates an excited state, even if one begins in the Bunch-Davies vacuum. We leave for the 
future a more systematic exploration of this connection. 

Finally, it would be most valuable to develop a broader understanding of the connection, 
if any, between symmetries and signatures in models of large-field inflation. 
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A Notation and Conventions 

In this appendix we review our conventions, emphasizing differences with the existing litera- 
ture. 

A good starting point is the ten-dimensional string- frame action 

5*10 — ^27r) 7 a; /4 j < ^ X ^/^ strin ^ ^ e ^ '^string ~ ^l^ 2 ^ (A-l) 
24 Rcmember that 2k\ = (2n) 7 a' 4 . 
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which after the rescaling to the ten-dimensional Einstein-frame metric e ^^ 2 g s tring,MN 
9e,mn becomes 



S 



1 



10 



(2vr)V 4 



d w Xy/gE [ Re - ^g s \dC 2 



(A.2) 



where we assumed that the axio-dilaton is r = i/g s - Upon compactifying on a six- dimensional 
manifold Y, the resulting four-dimensional reduced Planck mass is 



p (2vr) 6 a /3 q'tt a'rr ' 



(A.3) 



where Ve is the (dimensionless) Einstein volume of Y measured in units of 2n\/a'. When Y 
is (conformally equal to) a Calabi-Yau space, we have 



1 f J A J A J 1 f , ^ [ W; A Wj A L) K 1 , 7 K 

6 ( 27 r) 6 a'3 = ^ ^ ( 27 r) 6 a'3 = ^ ^ 



(A.4) 



where ouj for J = 1, . . . , ft, 1,1 are a basis of the cohomology H 2 (Y, Z) normalized such that 

ojj = {27c) 2 a'5j I (A.5) 



for a basis £/ of the dual homology H2(Y,/L). With the ansatz for the ten-dimensional RR 
two-form 



C 2 = —c(x)uj 

271 



(A.6) 



for some base two-cycle ui, we get a four- dimensional axion c(x) with periodicity 2n, as 
can be seen e.g. via S-duality starting from the world-sheet coupling J Bj(2-noi \ The axion 
decay constant of c is 



9s 



f UJ A *u 



„/3 



Ml 12V e (2tt) 2 [ (2vr) 6 a 
The four-dimensional M = 1 Kahler potential for the Kahler moduli is 

K = -2\ogV E . 



(A.7) 



(A.8) 



B Induced Shift of the Four-Cycle Volume 



In this appendix we address the issue raised in §6.2[ the inflationary energy can correct the 
warped volumes of four-cycles in the compact space, leading to corrections to the moduli 
potential, and hence inducing new terms in the inflaton potential itself. 



25 Note that this choice differs from that in [5], where the axion periodicity was (2ir) 
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More specifically, if an NS5-brane wraps some cycle S, then a nonvanishing integral 
J s C2 7^ leads to the presence of energy that is localized near £ in the compact space; 
this energy corresponds to the increased tension of the NS5-brane. Moreover, there is a corre- 
sponding induced D3-brane charge via the coupling J C2 A C4. The increased tension creates 
a backreaction on the metric (and in particular, on the warp factor) of the compact space, 
while the induced charge sources five-form flux. We must determine whether these effects sub- 
stantially correct the nonperturbative effects that are responsible, in our KKLT-like scenario, 
for stabilization of the Kahler moduli. 

Whether the nonperturbative superpotential arises from gaugino condensation on D7- 
branes or from Euclidean D3-branes, it is exponentially sensitive to the warped volume of the 
four-cycle wrapped by these D-branes. Therefore, we will carefully consider the possibility of 
an inflaton-dependent shift of the warped volume of various four-cycles. 

Concretely, we will consider a fivebrane/anti-fivebrane pair wrapping two homologous 
cycles, and will compute the leading correction to the volume of a particular four-cycle in the 
same throat region as the fivebranes. This will serve as a conservative upper bound on the 
effect of the worldvolume flux, as more distant four-cycles would be more weakly affected. 

It would be very interesting to perform a systematic study of this backreaction in the four- 
dimensional effective theory and in ten-dimensional supergravity/string theory. We leave this 
task for future investigation. In what follows, we simply show that the effect described above 
can be ameliorated by choosing an appropriate configuration. 

The are two mechanisms to suppress the backreaction on a given four-cycle volume. A 
first improvement comes from choosing a setup in which the leading backreaction is due to 
a dipole as opposed to a monopole potential. This allows for a parametric suppression. The 
second improvement can be achieved by a carefully chosen geometry of the four-cycle under 
consideration. In general, this latter mechanism requires fine tuning. 

The problem of estimating the backreaction from two-form flux on an NS5-brane pair 
may be simplified by a series of approximations. First, the inflaton-dependent backreaction 
is generated by the increased tension and the induced D3-brane charge of the NS5-branes, 
which may be understood as corresponding to some number of D3-branes (or anti-D3-branes) 
dissolved in the NS5-branes. In practice, it is much simpler to study the effect of the D3-branes 
themselves; this captures the leading inflaton-dependent contributions. 

The configuration of interest involves an NS5-brane wrapping S, with 

/c 2 = iV w , (B.l) 

as well as a distant NS5-brane wrapping a homologous cycle but with opposite orientation. 
(We will refer to the latter object as the anti-NS5-brane.) Next, we recall that the COBE 
normalization requires each fivebrane to be in a warped region. Let us denote by N/2 the 
amount of D3-brane charge that creates the background warping for each of the fivebranes,^ 

26 More generally, one could consider different degrees of warping for each fivebrane; extending our consid- 
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In light of the above discussion, we may approximate the fivebrane by a stack of N w 
D3-branes and the anti-fivebrane by a stack of N w anti-D3-branes. Combining this with the 
background D3-brane charge, we conclude that a convenient proxy for our system consists of 
two stacks of D3-branes, which we call A and B respectively. The first consists of N/2 + N w 
D3-branes and the second of N/2 D3-branes and N w anti-D3-branes, which we may more 
conveniently represent as N/2 — N w D3-branes and N w brane-antibrane pairs. 

Next, using the results of [H], we recognize that the leading backreaction effect comes from 
the total D3-brane charge on each stack, while the brane-antibrane pairs lead to subleading 
effects that are suppressed by powers of the warp factor. Thus, we can simplify even further, 
so that at last we are considering a supersymmetric system involving two stacks that contain 
N/2 + N w and N/2 — N w D3-branes, respectively. 

Equipped with this much simpler system, we may now estimate the inflaton-dependent 
backreaction, by computing how the presence of the stacks A and B leads to a A^-dependent 
change in the warped volume of some four-cycle. 

We choose the usual D3-brane ansatz 

dsl = JW=Hy)dsl+ y/H^jdsl (B.2) 
p 5 = (1 + *)dH^ x A Vol 4 , $ = const. (B.3) 

The resulting equation of motion is linear in H(y). Therefore we may simply add the solutions 
obtained in the presence of either of the two individual stacks. Once the resulting warp factor 
is used to compute the volume of a four-cycle, the N w dependence of the result gives us an 
estimate of the inflaton-dependence of the nonperturbative superpotential. 

We tackle the problem in two steps of increasing complexity. First, in §B.l| we give a very 
simple, (conformally) flat toy example in which the calculations are easy. This already shows 
the relevant features of the more complicated solution. The inflaton-dependent shift of the 
volume can be suppressed by having the distance between A and B much smaller than the 
distance between the four-cycle and either of A and B. This corresponds to a configuration 
in which the leading interaction is via a dipole. In addition, one can fine-tune the four-cycle 
embedding so that the iV^-dependent correction to its volume actually cancels, 
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Then, in £B.2 £B.2.1 and §B.3 we describe the case of a resolved conifold using the 
solution of [521 E3]- We consider a particular holomorphic embedding of a four-cycle and 
compute numerically the inflaton-dependence of its warped volume. 



B.l A simple illustration of the suppression mechanism 

Consider two stacks of N/2 ± N w D3-branes, called A and B, respectively, in conformally flat 
space M4 x IR 6 . The A stack is located at the origin of IR 6 and the B stack is located at some 

erations to this case is straightforward. 

27 Although suppression from symmetry of the embedding appears unappealing because of the fine-tuning 
required, one should keep in mind that it could conceivably be enforced by a discrete symmetry of the 
compactification . 
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Figure 9: This diagram illustrates the positions of the A and B stacks of D3-branes in M 6 , 
the choice of the angular coordinate 9, and, in blue, a (topologically trivial) four-cycle. 



position (u, 0, 0, 0, 0, 0) foi 28 u G M + , where we have chosen spherical coordinates (see figure 



B.l) with the metric 



dsl = dr 2 + r 2 (d6 2 + sin 2 9dQl) 



(B.4) 



where dfl 4 is the volume form of S A . With the usual D3-brane ansatz (B.2) one finds the 
solution 



H 



it a rip, 
1 H — H , 

r 4 l r 2 _|_ u 2 _ 2 ru cos Q\2 ' 



R 



r 

4:7rg s a' 



Let us consider a (topologically trivial) four-cycle S 4 defined b> 29 
unwarped volume is V 4 



H(r, 9) sin 4 6dVL 4 



2|-r sin 9. The warped volume is 



(B.5) 
(B.6) 

fi and 9 = 9, whose 
(B.7) 



^4 



V 4 



1 + — ll + O - 

I 1 \ \P 



11 f -11 I 

N, n — I 4cos9 — 2 — yO 



28 We choose the letter u in analogy with the setup of the next subsections, where the distance between the 
A and the B stack is given by the resolution parameter of the resolved conifold. 

29 We choose the letter fi in analogy with the (usually complex) parameter appearing in other known 
embeddings of four-cycles in the conifold [5H |5S] . 
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From this result, one can see that a fine-tuning of the embedding can suppress the backre- 
action, i.e. if cos(9 ~ u/(2fi). On the other hand, a parametric suppression is also clearly 
visible. Both the factors N w /N and u/fi can be made small by construction. The former 
must be small in order for the background geometry to be at all trustworthy. The latter can 
be made small by arranging for all the four-cycles bearing nonperturbative effects to be far 
away in units of the separation of the fivebranes. Physically, this means that the four-cycle 
is sensitive only to the dipole field generated by the A and B stacks. 



B.2 The conifold and its resolution 



In the following we review some relevant definitions and conventions regarding the resolved 
conifold. The treatment is based on [56, 53J. The (singular) conifold is a cone over T 1,1 
(the coset space SU{2) x SU(2)/U(1), which is topologically S 2 x S 3 ). It is defined as the 
hyperspace in C 4 that is a solution of the complex constraint 



detW = det 



X U 
V Y 



XY - VU = 



(B.8) 



where (X, U, V, Y) are coordinates on C 4 . The resolved conifold can be defined as the zero 
locus in C 4 x CPx of the two linear complex equations 



X U\ X t 
V Y H A 2 







(B.9) 



where (Ai, A 2 ) are complex coordinates on CPi, i.e. they are identified by (A]_, A 2 ) — (otAj, aA 2 ), 
for every a 6 C». For every W ^ (B.8) and (B.9) are equivalent, but when W = 0, i.e. at 
the tip, (Ai,A 2 ) are arbitrary and (B.9) defines a CPi ~ S 2 . The radial direction is defined 
by 



(B.10) 



One can check that the resolved conifold is an 0(— 1) © 0(— 1) bundle over CPi with fiber 
C 2 . If we define A = A 2 /Ai, then we can choose coordinates on a patch H + = {A ^ 0} of the 



resolved conifold using the following solution of (B.9) 



W 



-XU U 
-XY Y 



(B.11) 



Defining A = Ai/A 2 , one can find coordinates on a complementary patch H 
the following solution: 



W 



X 
V 



-XX 
-XV 



{X 7^ 0} using 



(B.12) 



54 



The complex structure is given by 

= dUAdY AdX = dV AdX AdX. (B.13) 

For later use, we introduce a parametrization of the resolved conifold in terms of real coor- 
dinates and give the explicit Kahler metric. We start by noting that a particular solution of 
(B.9) is given by 

r \ I \ 



w °={ oj- { o J =* A = ' < R14) 

The base of the resolved conifold with respect to r can be obtained by acting on this solution 
with two SU(2) transformations, Li and L 2 , 



cos $i e i(*+*) _si n | e -|(^-^) 
sinfei^-^ cos|e-5 { ^+^ ) 



i = l,2 (B.15) 



written in terms of Euler angles. This gives 



W = L 1 W Ll (J* j=L 2 ( A q ° J (B.16) 

which depends only on the combination ip = ipi + ip2- A Kahler metric on the resolved conifold 
with resolution parameter u is given by [53] 

ds\ = K-\p)dp 2 + l -^p)p 2 el + \p\e\ + eJJ + \{p 2 + Gu 2 ) (e 2 , + ejj , (B.17) 

where 

o 2 + 9m 2 

«M = ^- (B-18) 
Here, following [53] . we have defined a new radial coordinate p by 

r 4 4/fV + 6« 2 ) , (B.19) 



9' V3 

The explicit expression for the e's is 



ef = dip + } j cos 6jd(j)i , e 6i = d9i , = smOidfa . (B.20) 
i=i 
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B.2.1 The XUY embedding 

In this subsection, we consider a particular holomorphic embedding of a four-cycle. A simple 
embedding would be A = p because this is trivial to solve for in real coordinates, tan 6 2 = p 
and (f>2 = for /jGl The trouble is that this embedding reaches the tip, and in fact r is 
unconstrained. This can also be seen from 

r 2 = (1 + |A| 2 ) (\U\ 2 + \Y\ 2 ) . (B.21) 

As a result, this embedding does not give us the dipole suppression factor analogous to the 



(u/ n) of appendix B.l The next-simplest embedding (whose defining equation depends on 
r) is 

\UY = p 3 , fieU, (B.22) 

which in real coordinates gives 

3 3 

V> = 0, sin(0 2 )sin(0i) = 4^- ~ ^ for large r . (B.23) 

r 2 p A 

After some algebra (in particular, expressing d6 2 as a function of d6\ and dr) one finds 
the metric in terms of dp, d9 2 , d(f>i and d(j) 2 . Its determinant g™ d is independent of </>i j2 and 
reads 

ind = p 2 csc 4 (fl 2 ) x 

20736 (9m 2 + p 2 f (p 4 (9w 2 + p 2 ) sin 2 (<9 2 ) - 54/i 6 ) 

(6m 2 + p 2 ) (9u 2 p + p 3 ) 2 cos (40 2 ) 

-4 cos (20 2 ) (486w 6 p 2 + 189w 4 p 4 + u 2 (24p 6 - 324p 6 ) - 27p 6 p 2 + P § ) 

+3 (54p 2 (9m 6 + 2p 6 ) + 189w 4 p 4 + 864w 2 p 6 + 24w 2 p 6 + p 8 ) ) 2 . (B.24) 

We see that there is a boundary beyond which the sign of the determinant becomes 
negative, which thus defines the integration boundary in p, # 2 -space: 



Pmin(0 2 ) = V3p \lA-j^- (i-^) ^ ^3-21/3^ csc 1 / 3 ^) for ^ < 1 



with : A 



a y - ^ 2 ; - uea, y - - 



'-csc 2 (# 2 ) ( 2^-csc 2 (0 2 ) J -^ + csc 2 (0 2 ) . (B.25) 



B.3 The shift of the four-cycle volume 

The solution with the branes smeared over the S 2 was obtained in [53]. Later, the solutions 
with pointlike sources were given in [52] . If the D3-brane stacks are at the north and south 
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pole of the resolution S 2 , respectively, i.e. 9 2 = vr — 9 2 = 0, then one finds 

H = ^(2Z + l)i^(p) [^(cos(^)) + L^(cos(0 2 ))(-l) ? ] 



H, 

La,b 



2 

9m 2 p 2+ 2/ 3 

(3n) 2 ^r(l + /3) 2 



2 Fi A 1 + /?, 1 + 2/5, 



9m 2 



r(i + 2/3) 



/3 = v 7 ! + (3/2)/(Z + 1) 
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(B.26) 

(B.27) 
(B.28) 
(B.29) 



where 2 Fl is a hypergeometric function. We want to integrate this warp factor on some 
supersymmetric four-cycle £4. This gives us an estimate of the inflaton dependence of the 
gauge kinetic function of a stack of D7-branes wrapping £4. 



Using this information and (B.26) and (B.27) we can now calculate the integral 



V, 



warped 



dp d9 2 



£ 4 



g? d H(p,9 2 )=4n 2 dpd9 2 J-g? d H(p,9 2 ) (B.30) 



numerically, as a function of p. To facilitate this we will expand (B.26) up to 
dipole term, and take the large p limit 



H(p,9 2 ) 



V 



N^v 2 
N p 2 



1 + 3(2^ + l)4^P € (cos0 2 ) 



H i=0 (p) + 6H e=1 (p,9 2 



We can now calculate 



(o) 

warped 



^Lalped\ uu 2l 



ft fPR 

4tt 2 / d9 2 / dp\/-gi nd H £=0 (p) 
An 2 d9 2 dp\/-gl nd 5H e=1 (p,9 2 + 59 2 ) 

JO J pminifil) 



1, the 
(B.31) 

(B.32) 
(B.33) 



where pr 1 denotes a UV cutoff to compactify the resolved conifold geometry for the 
purpose of integration, and 59 2 denotes the angular misalignment of the D3-brane dipole 
configuration with respect to the four-cycle symmetry axis at 9 2 — 7r/2. 

As \J #i nd is a symmetric function with respect to 9 2 = n/2 and 5He = i(p,9 2 ) is anti- 
symmetric with respect to 9 2 = n 2 , we immediately find 5V war ped(8@2 = 0) = 0. So by 
fine-tuning a Z^-symmetric configuration we can forbid the i = 1 term in the warped volume, 
whose corrections in this case start with the £ = 2 quadrupole terms. 

We will now display the numerical results for the case S9 2 = —ir/4 in which the i = 1 term 
will not vanish under the integral, and compare the scaling with p between SV warpe d(S9 2 ) and 



V wiped- This is displayed in Fig. 



- const, while the 



10 



We see clearly that the leading £ = term scales as 
1 dipole term scales as (u/p) 2 . Therefore, the i = 1 dipole term 
has a parametric suppression {u/p) 2 relative to the leading I = term, and can therefore 
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Figure 10: 1st row: Plot of d and 5V warpe d((56 l 2) as functions of p at constant w = 0.01. 

2nd row: Plot of V^ rped and 5V wa r P ed(SQ2) as functions of w at constant p = 0.1. The leading 
£ = term scales as = const, while the £ = 1 dipole term scales as (u/p) 2 . Note that 
the I = scaling ensues only in the strictly noncompact limit (i.e. when the integration goes 
all the way p — > oo), while for a finite cutoff, resembling a crude approximation to a compact 
setting, there remains a weak dependence of the £ = term on /j, of the form (u/p) s , where 
5 — > for pfe u ifc — » oo. For the example we have chosen e = N w /N = 0.1. 



be made parametrically small (even in the non-Z2-symmetric general situation) in the limit 
where the four-cycle recedes far from the resolution S 2 (i.e. in the limit of large fi/u). 

Let us finally note that this relative suppression of the I = 1 term with (u/fj,) 2 might have 
been guessed without any integration, as the integration boundary tells us that p m j n (6 , 2 ) > 
Pmm(7r/2), which corresponds to r > or p > fj,, and thus the relative scaling u 2 / p 2 

should be replaced by the scaling u 2 /p 2 . 

C The Kaluza-Klein Spectrum 

In this appendix we obtain the (5+l)-dimensional effective action for a D5-brane wrapped on 
a two-cycle with J B ^ 0. We show how a Kaluza-Klein reduction to four dimensions leads 
to masses that are suppressed with respect to the fluxless case. We then comment on the 
consequences of these light KK modes for axion monodromy inflation. 
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C.l The effective theory 

The DBI action for a D5-brane is 

S = T 5 J d A x dy dz^-det (Gff + F ab ) , (C.l) 

where y, z are two coordinates in the internal space, which we take to be toroidal for the 
purpose of this derivation. The indices are defined as follows: worldvolume indices are a,b — 
0, . . . , 5; spacetime indices are /j,, v — 0, . . . , 3 as usual; ten-dimensional indices are M, N = 

0. . . . , 9; six-dimensional compact indices are m, n = 4, . . . , 9; and indices transverse to the 
D5-brane are i,j — 6, . . . , 9. We first expand the square root using 

v / det(M + 5M) = y r ^mAl + -i:i{M Q 1 5M) + -[Ti{M ( ; l 5M)} 2 (C.2) 

[2 8 

-^{M^5MM^5M) + . . . | . (C.3) 
We will consider a background with two-form flux on the two-cycle 

J T= J B = J dyA dzB yz (x, y, z) = b(x) = b , (C.4) 

1. e. the four-dimensional axion field b(x) has a homogeneous expectation value that is ap- 
proximately constant, up to terms suppressed by the slow-roll parameters. So the background 
is given by 

BmN = bSMy^Nz — b5Mz$Ny , Fab = , (C5) 

ds io = Siivdx^v + g yy dy 2 + g zz dz 2 + 2g yz dydz + gijdy l dy 3 . (C.6) 



Hence 



(^A»U= j g yy g yz + b ) . (C.7) 

g zy -b g z 



izz 



The perturbations are 

(5M) ab = d a X%X^g tJ + B tj ) + F ab + 5B ab . (C.8) 

The calculation is simplified by the block-diagonal form of the background M . The 2x2 
block is the sum of a symmetric and an antisymmetric piece that we call S and A respectively. 
We have that 

det(A + S) = det(A) + det(S) , (C.9) 
{b+A} - b det(A) + det(S) + A det(A) + det(S) ' 
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which substantially simplifies the calculation. Using (C.2) we get at leading order 



S 



T 5 J d 4 x dy dzy/-g 4 \/g 2 + b 2 

1 92 



1 + -d^d^Xi 



+ 



2g 2 + b 2 
1 2b 
'2g 2 + b 2 



(dyX'&Xi + dzXWXi) 
(d y X i d z X i SB ij + F yz + 5B yz ) 



(C.11) 
(C.12) 
(C.13) 



g£„. After a KK reduction one finds the four-dimensional kinetic and 



where g 2 = g yy g Z ; 

potential terms, in the first line, as well as the Kaluza-Klein mass terms, in the second line. 
The Kaluza-Klein masses in the presence of fluxes are 



m 



92 



bKK 



92 + b 2 



m 



KK > 



(C.14) 



where ttikk are the Kaluza-Klein masses in the absence of fluxes. This leads to the central 
point of this appendix: for b ^> 1, the Kaluza-Klein masses are suppressed by a factor of 
> ~ L 2 , 



< 1 



30 



This phenomenon is intuitively understood in the T-dual picture in 
which flux becomes the angle of the D-brane. A large flux means that the T-dual brane winds 
around the torus many times, and thus becomes quite long. The Kaluza-Klein reduction 
of the fields living on the worldvolume of the T-dual brane therefore produces b-suppressed 
Kaluza-Klein masses. 



C.2 Effects of the light Kaluza-Klein modes 

Throughout this paper we have been careful to work in parameter ranges for which the typical 
Kaluza-Klein mass scale tukk obeys 3> H, as required for a consistent four- dimensional 



analysis of inflation. However, from (C.14) we learn that a subclass of Kaluza-Klein modes, 
namely those associated with transverse excitations of the fivebrane, have considerably smaller 
masses, tti^kk *C itlkk- For the numerical examples we have considered, we find that, very 
roughly, tti^kk ~ (fc/ f)H, where f c is a fiducial value of the decay constant, f c ~ 10~ 2 M p . 
Therefore, for constructions with small values of /, the transverse excitations of the fivebrane 
can be lighter than H. 

We leave a comprehensive study of this constraint for future work, as a proper implemen- 
tation plausibly requires a more explicit compact model that we have been able to present in 
this work. In particular, one should carefully compute the Kaluza-Klein mass, incorporating 
anisotropy in the geometry, warping, and, as we have explained above, the effect of world- 
volume two-forms. To accomplish this, one needs a reasonably explicit construction of the 
warped throat region, of the two-cycle within the throat, and of the gluing of the throat into 
the compact space, which are beyond the scope of this work. 

30 We have assumed for simplicity that the internal space is isotropic, with typical size L\fct! '. 
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In this appendix, we will restrict ourselves to some qualitative statements that explain 
how our inflationary analysis can be consistent even in parameter regimes for which m^xx is 
slightly smaller than H . Broadly speaking, one might worry about corrections to the inflation- 
ary Lagrangian, and about new contributions to the cosmological perturbations. Concerning 
the first point, we remark that the excitations of the fivebrane depend on the inflaton expec- 
tation value only through their masses. Therefore, the primary correction to the background 
evolution from these light modes would come if large numbers of Kaluza-Klein particles were 
produced by the time-dependent background. In practice, the particle production is negligi- 
ble, as can be seen by computing the adiabatic parameter fhbKK / m< t>KK an d substituting the 
constraints on the volume, and hence on the Kaluza-Klein mass, from £j6j 

More generally, let us stress that only a small subclass of the Kaluza-Klein modes (a 
small portion of the tower of excitations of the fivebrane) have masses smaller than H. From 
the viewpoint of the inflationary analysis, these fields constitute a small number of harmless 
spectators. These light fields will fluctuate, absorbing energy, but this yields a very small 
correction unless the number of fields approaches (M p /H) 2 . Moreover, any entropy pertur- 
bations produced by these fields can turn into visible isocurvature perturbations only if their 
decays are distinct from that of the inflaton. Although we have not specified a concrete 
reheating mechanism, one can argue that the most straightforward scenario involves visible 
sector degrees of freedom that are well-separated in the compact space from the inflationary 
fivebranes. Thus, we expect that excitations of the fivebranes will not give visible isocurva- 
ture perturbations, because they must first decay [53 EE] to degrees of freedom localized in 
the inflationary throat, just as the inflaton does, and will plausibly do so with rather similar 
couplings, as the modes correspond to small excitations of the NS5-brane that drives inflation. 



D Numerical Examples 

In this appendix, we specify two different sets of intersection numbers and show the relevant 
formulas for the volumes. For these two toy models, we explicitly performed the moduli 
stabilization outlined in §6.5[ finding numerical values leading to the dot in figure [8] 



D.l Intersection numbers: set I 

We consider as a toy-model Calabi-Yau manifold one with = span(tu i , u + ) for the 

orient if old-even homology two-cycles and Hj 1 = span(c<j~) for the orientifold-odd homology 
two-cycles. We assume the following simple set of intersection numbers 

clll = c LL+ = c + __ = 1 , (D.l) 

with all the others vanishing. We believe that, although very simplistic, the above toy model 
captures the relevant features of more realistic constructions. Notice that the intersection 
numbers in a basis for the homology of the covering space of the orientifold, i.e. without a 
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definite parity with respect to the orientifold projection, are just linear combinations of those 
given above. 

Using the standard relations 

V E = ^c a/3 yV a v p v' r , r a = d v oV E = -Go/^iA; 7 , (D.2) 



one finds 



and 



v L = V2T+, v + = , (D.3) 



V E =^r L -^-rl'\ (D.4) 
D.2 Intersection numbers: set II 

Again assuming H+ 1 = span(c<j L , u + ) and ifi' 1 = span(c<j~), we consider the intersection 
numbers 

clll = c L++ = c + __ = 1 , (D.5) 
with all the others vanishing. We find 

*L = -j= ((TL + T+) 1 / 2 + (r L - T+ )W) = ^ ((T L + T + )W ~ (t L - T + )^) , (D.6) 



and 



V E = ^((r L + r + )^ + (r L -r + )^) . (D.7) 
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